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Preface 



This thesis is a compendium of our work on extension of basic Skyrme model to 
arbitrary representations of SU(2) group, hoping that higher representations would 
be helpful for more adequate description of static baryon properties jl], ^ |^] . 
General ideas and historical remarks. The idea that the ordinary proton and 
neutron might be solitons^ in nonlinear model has a long history. The first sug- 
gestion was made by T.H.R. Skyrme about 40 years ago J5J. The essential feature 
of the theory is the representation of the fundamental field quantities in terms of 
angular variables rather than linear ones. Realistic three-dimensional model is pos- 
sible only when there are also three angular variables. The condition is satisfied by 
the pion fields of nature. The periodicity of angular variables introduces a new con- 
stant of motion, which measures the number of times that space (three dimensions) 
is mapped by the fields onto the elementary volume of angular space and which 
can be interpreted as a baryon number. The origin of the new constant of motion 
is related to topological features of the Skyrme model. In contrast, conservation of 
energy and momentum follows from space-time symmetry, as usually. 

The mathematical construction outlined above is to ensure possibility of particle- 
like states "of a kind that cannot be reached by perturbation theory and which 
cannot necessarily be discounted by general arguments" Q|. To provide readers a 
link between fundamental theory of strong interactions (QCD) and Skyrme model 
we need to consider briefly the chiral symmetry concept and, therefore, the idea of 
isospin. 

The concept of isospin was introduced in 1933 by W. Heisenberg, who consid- 
ered proton and neutron as different projections of single state^. From contem- 
porary point of view W. Heisenberg actually assumed SU(2) (flavour) symmetry 
of nuclear interactions. In 1962 M. Gell-Mann succeeded much more in suggest- 
ing very predictive SU(3) (flavour) symmetry of strong interactions (The Eightfold 
Way) and the concept of isospin was extended to all baryons. The SU(3) symmetry 
had enormous influence in becoming of QCD: it was realized that each basis ele- 
ment, i.e. a product of quark and anti-quark functions, can be identified with some 
hadron state. The entire basis, therefore, is interpreted as multiplet of hadrons, 
belonging to irreducible representation of (flavour) SU(3) group. 

A revival of interest in the Skyrme model [jS], |<|, begins from the work Q of 
G.S. Adkins et al., who demonstrated that this model could fit observed properties 



2 Soliton history begins in 1834, from D.S. Rassel's (1808-1882) "great solitary wave". There 
have been, however, no more than twenty scientific works during the period 1845-1965, directly 
related to solitons |4|]. 

3 W. Heisenberg even suggested to explain interaction between proton and neutron by particle 
exchange Q. The existence of pion, however, was predicted by Yukawa theory in 1934. The 
particle was discovered by G. Lattes, H. Muirhead, G. Occhialini and S.F. Powell in 1947. 
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of the baryons to an accuracy of about 30%. This rebirth of attention was stimu- 
lated by the belief that some such model is the long-wavelength limit of QCD, as 
reviewed, e.g., in Refs. [iofl . The interplay between various phenomenological 
models and QCD is still open problem 



11 



Extensions and modifications of basic Skyrme model. In the intervening 
period, there has been a large number of works extending the range of applica- 
tions, modifying and extending the model, and improving the way in which conse- 
quences are drawn from it. This work also serves as an extension of basic Skyrme 
model to arbitrary SU(2) representation. Among the further applications, the most 
prominent have been to pion-nucleon scattering [12, 13 1 and the two-nucleon prob- 
lem [14, 15|, and, very recently, to multi-soliton "chemistry" [16, 17, 18]. The 
basic Skyrme model has been extended in various directions. Excluding the men- 
tion of models that contain the quarks explicitly, one encounters in the literature 
models with higher-order terms involving the same fields [19, 20 1 or even higher 
unitary groups [21], models in curved space [22, ^| and models in which vector 



25[ |26[ , as well as extensions to include strange |2J 
Related but more fundamental extension of the 



mesons have been added |]24|, 
and even charmed mesons 
basic Skyrme model is the incorporation of the Wess-Zumino term into this the- 
ory. This term eliminates an extra discrete symmetry that is not a symmetry of 
QCD [29] and, therefore, has far reaching consequences and, of course, no lack of 
attention |30|. 

All these models are first presented as classical field theories, since one can 
do much physics using only selected classical solutions. The need to address the 
problem of quantization is, however, manifested in the intrinsic properties of the 
classical solution. 

The capability of extracting interesting physics from the Skyrme model is 
grounded on the existence of a special solution of the classical field theory, the hedge- 
hog skyrmion. Like all interesting classical (or mean field) solutions it breaks some 
of the symmetries of the underlying Lagrangian. The hedgehog skyrmion violates 
translation, spatial rotation, and iso-spatial rotation symmetry. The restoration 
of these symmetries requires, at the very least, the quantization of the generators 
of the symmetry transformations and of the associated canonically conjugate col- 
lective coordinates. As a consequence, maximum attention has been paid to this 
aspect of the problem of quantization [31 1. In addition, to study pion-baryon scat- 
tering 1 12, 32, 33 L it is necessary to discuss quantization of the small oscillations 
of the pion field [34] (for different approach see ]35|). There have also been some 
discussion on quantization of radial oscillations 36, 37 , 38] in connection with 
problems of stability. Other quantization methods applied to Skyrme model in- 
clude cutoff quantization |39|, (which uses short-distance cutoff e : F(e) = n), 
the general covariant Hamiltonian method [40, 41, |42] | (which preserves the orig- 
inal symmetry of classical Lagrangian) , Kerman-Klein quantization procedure |43| 
(based on formal quantization of entire classical field). Due to rich and beautiful 
mathematical structure the model has numerous applications. 
Applications of Skyrme model. Despite the original model has been introduced 
to describe strongly interacting particles, there are attempts to apply similar gauged 
construction to describe weak interactions ("electroweak skyrmions") [44| (p. 250). 
Apart from high energy physics the model proved to be useful in cosmology |45| 
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and solid state physics^] ||46| , 47]. Before brief review of the manuscript organization 
we explicitly state main tasks this thesis is intended to solve. 
Main tasks: 

1. Investigate representation dependence of the quantization procedure |48| 
applied to the SU(2) Skyrme model. 

2. Numerically evaluate obtained expressions for physical quantities and com- 
pare the results with experimental data. 

Scientific novelty. This work demonstrates the new possibilities to extend basic 
Skyrme model to arbitrary representations. Quantization of the Skyrme model (in 
collective coordinate approach from the outset) yields different quantum Lagrangian 
density for each SU(2) group representation j. The classical limit of these quantum 
Lagrangian densities is the same original Skyrme Lagrangian density. For the first 
time it has been shown that stable quantum solitons exist both for spin, isospin 
I = 2 and I = | states. These quantum solitons possesses Yukawa asymptotic and, 
therefore, imply non-vanishing pion mass. Noether currents, magnetic momenta 
etc., operators have been calculated and numerically evaluated in this approach for 
self-consistent quantum chiral angles in various representations j. 

The generalization considered in the work has far-reaching consequences and, 
we believe, can readily be extended to other models and theories. 
Manuscript organization. The manuscript is organized into three chapters plus 
appendices, containing numerous tables and illustrations. Chapters and sections 
(if structure of the latter is complicated enough) have short information about its 
content and, therefore, not need to be repeated here. We find useful, however, 
briefly to describe what purposes each chapter is intended to serve. Chapter I 
contains mathematical formulation and physical motivation of the Skyrme model 
in background level. Apart from few presentation details it contains no new re- 
sults. We give formulation of classical Skyrme model in group theory terms and 
introduce mathematical apparatus which is convenient for model formulation in 
arbitrary representation in Chapter This chapter includes results of Ref. |49|. 



Chapter fl| contains main new results and deals with the quantization of the 
Skyrme model. 

Despite the thesis has no lack of references when investigating concrete prob- 
lems we found useful to provide a list of sources about the entire model. These are 



books |44|, 50| and review articles |5lj, |52|, |53|]. Literature on solitons currently is 



untraceable^, but we still mention few books, namely 0, |55) to begin with 
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CHAPTER I 



Introduction to the Skyrme model 



This chapter is intended to provide very short but more or less consistent 
introduction to the Skyrme model. From this point of view it is essential to clear out 
difference between models realized linearly and nonlinearly. The simplest examples 
are linear and nonlinear a models. Skyrme model then arises naturally by adding 
the fourth-order term in field functions to nonlinear a model Lagrangian. This 
term enables existence of stable soliton in three spatial dimensions (skyrmion) and, 
therefore, is called a stabilizing term. 

1. Linear a model in two spatial dimensions 

We start with linear a model. When physical boundary conditions are imposed, 
all model solutions fall into disconnected classes regardless of what equations of 
motion are. It is this property which is peculiar to nonlinear models only and play 
an important role in the Skyrme model particularly. 

1.1. The Lagrangian. Let's take $ to be a scalar doublet of real fields 
(<I>i,$2) and consider the Lagrangian 

JSf = -ld M $ r d"$ r -(|*| 2 -a 2 ) 2 ; r = l,2; a 2 > 0. (I.l.l) 

Here we understand |<fr| 2 = <& r $ r and \i — 0,1,2 (0 denotes time component). 
Assume that for x — ► oo field configuration tends to some constant statefj] 

$ 2 + $ 2 ^a 2 , for|x|->oo (1.1.2) 

with approach rate, which guarantees finiteness of total system energy E. Then the 
set 4>oo = ($001, $002) of all fields at spatial infinity make up a circle S . Spatial 
infinity in argument plane also can be imaged as a circle S , with infinitely large 
radius 

x\+xl=R 2 , R^oo. (1.1.3) 

Field $00 maps a circle S 1 to a circle S 1 , $00 : S 1 — > S 1 . Identification of 
infinities with a circle of infinite radius does not involve any topology change. The 
hint is similar to coordinate system change. We have here spaces — both argument 
and function space — flat. Moreover, these spaces are vector spaces as well. Con- 
sider first a vacuum (or trivial) solution #vac(^i> ^2) = (a,0). At spatial infinity 
the field (<&oo)vac = lim <& vac maps all points of S* 1 circle (argument infinity) to 

R — >oo 

the same vacuum point (a, 0) of function space S 1 . Thus (<froo)vac is characterized 
by zero winding number: (3?oc)vac € Qo°- To (3>oo)vac we can, therefore, associate 

1 Only these states are interesting from physical point of view: energy at infinity should be 

zero. 
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Figure 1.1. Example of zero winding number maps. 

all maps which are homotopic to (<froo)vac- The set of all maps makes up 
the trivial sector Qg° of system configuration space Q°° (see Fig. [Tl| ). Of course, 
one can choose any other vacuum state simply redefining <i? vac (x) = (a/v2, a/^/2) 
or more generally 3?vac(#i, £2) = (acosx, asinx), where x is fixed from [0,27r] 
(see Fig. C ). All vacuuma belong to trivial n = sector and are homotopic to 
4* vac = (a, 0) . The homotopy can be defined as /12 (t) = (a cos tx, a sin t%) t S [0,1] 
(see Fig. [Oj^c). One could think about gauge freedom corresponding to transforma- 
tion from one vacuum to another. With any vacuum choice model SO (2) symmetry 
becomes spontaneously broken. Physical motivation comes from tunnelling possi- 
bility: vacuuma aren't separated by any potential barrier, therefore, in the case of 
such vacuuma interference fields must rotate everywhere in space. This involves 



infinite rotational energy [44] (p. 107). Of course such interference also restores 



SO(2) symmetry of the ground state and, therefore, should be forbidden. 

1.2. Soliton sectors and invariants of linear a mo del. A field (0) of 
the sector Qf° can be defined as &)^(8) — (a cos 9, a sin 9). When 9 runs from to 
2tt all points on S 1 are covered once and only once. is a typical winding number 
one map (see Fig. L2). The equivalence class of maps homotopic to forms the 



winding number one sector Qf. The winding number n map can be defined as 
^oa(@) = (a cos n9, a sin n9), where n is an integer, due to our requirement of 
* single valuedness on S 1 : #^ } (0) = ^(2tt). Q { n oo) consists of all *£ } 
homotopic maps. We remind that with any n (not only n = 0) satisfy (1.1.2). 



From this picture it becomes clear that it is not possible to deform a field 
to a field continuously, if n 7^ m. Indeed, for the action we need to cut the 

mapping curve and take off \n — m\ twists. Thus the fields in sector q! 00 " 1 are not 
homotopic to fields in <2m°^ f° r n ^ m. As a consequence, system configuration 
space Q°° falls into an infinite number of disconnected components = {^x}, 
Q°° being a union U n Q™. The same is then true for the space Q of fields defined 
over all space: Q = U„ Q n . Here Q n is the space of all configurations $ < - n - ) whose 
limit as |x| — * 00 is an element of . For example, $ < -™- ) at given time to 
can be defined by 

*(">(x,t )=/(|x|)*W(_l), (LL4) 

where /(|x|) is any smooth function, such that /(oo) = 1, /(0) = 0. 

The physical significance of the integer n associated with the field <fr^(x, t) is 
that it is an integral of motion. The integer is just the label of homotopy classes of 
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Figure 1.2. Example of winding number one map in linear a model. 



the fields at a fixed time. Sectors Q n (n ^ 0) are topologically stable in the sense 
that a field in Q n will not evolve in time to the vacuum (an element of Qq) or to 
a field in any other sector Q m (m ^ n), since evolution in time is a continuous 
deformation. 

The existence of inequivalent topological sectors leads to additional invariants 
in the theory. These new invariants of quite different origin is the most interesting 
and important point in such models. As a consequence, we have two kinds of 
invariants: 

• Invariants which are closely related to the symmetry of the system, under si- 
multaneous coordinate frame and fields change (corresponding to this frame 
transformations). We can find all these invariants by Noether theorem. 
Examples of the invariants are: energy, momentum, angular momentum, 
electromagnetic charge.^] 

• Invariants, involving boundary conditions in one or another way. Conser- 
vation of the number of particles in classical mechanics (in conservative 
systems) is an example. 

Let us concentrate on the second type of invariants. It has been shown that so- 
lutions from one sector cannot evolve in time to the solutions of any other sector. 
Consequently, we need discrete quantity to label each sector. The most natural 
choice seems to be the number of twists, describing mapping of one circle onto an- 
other. The number is called a topological charge. One could introduce a conserved 
topological current density, corresponding to the charge 

4 op ~ e^ p e rs d°$ r d p $ s ^,^ = 0,1,2; r,s = 1,2. (1.1.5) 



It is easy to check that expression dl.l.5|) has a divergence zero and meets our 



requirements. We stress that the current density is conserved irrespective of what 



2 The Lagrangian ( |[.l.l| ) is symmetric under rotations of SO(2). The SO(2) is known to be iso- 
morphic to U(l), therefore, theLagrangian Jz? = —^d^Qd^Q* -(<!>$* — a 2 ) 2 ; /i = 0,1,2; a 2 > 
can be chosen instead of ( |[.l.l| ). Conserved Noether currents exist corresponding to continuous 
symmetries (SO(2) or U(l)) of these Lagrangians. In the case of one complex field the conserved 
electromagnetic current has a simple form J° ~ i (<b*dn. & — (9 M $*<3?). These currents have noth- 
ing to do with the conserved topological current ( f .l.aj ). Also there is one interesting difference 
between two real and one complex field case. Namely, there is no real vector, which is invariant 
under SO(2) rotation, but there is a pair of complex eigenvectors (1 ± i) with eigenvalues e ±l ^ in 
complex plane (group U(l)). 
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the equations of motion are (because of the antisymmetric properties of e MCTp ) and 
thus it is topological. The corresponding charge density is 

J* op = divP, where P k ~ e k ie rs ® r di$ s ; fc, Z = 1,2. (1.1.6) 

The charge Q top ~ J jQ° p d 2 x is nonzero only for fields <fr with non-vanishing as- 
ymptotic^] |56| |, which is realized, for example, by the Higgs mechanism. 

1.3. Derrick theorem. The existence of topologically stable sectors and con- 
servation of topological current are independent on the Lagrangian form and, thus, 
on equations of motion. The presence of such sectors, therefore, does not guarantee 
that equations of motion actually have solutions in each sector. It is known that 



Lagrangian (1. 1.1) does not lead to nontrivial stable static solutions of equations of 



motion if only the spatial dimension is D ^ 1. This can be shown by simple scaling 



argument of Derrick [44, 57]. Suppose that $ c i is static solution and the energy 



= or (2-D)E 1 =DE 2 . (1.1.9) 

A=l 



of the solution consists of terms E — E\ + E 2 : 

Ei = \ J d B x(^(* cl ) s ) 2 and E 2 = y"d D xC/[((* cl ) s ) 2 ] i = 1, 2, . . . , D. 

(1.1.7) 

Under a scaling transformation 3? c i(x, t) — > $ c i(Ax, t) these terms scale as 

E = E(l) -> E(X) = X 2 -°E 1 + \- D E 2 . (1.1.8) 

Requiring that A = 1 corresponds to energy E minimum yields the condition 
dff(A) 
dA 

Since Ei,E 2 > 0, it follows that E\ = E 2 = when D > 2. This implies that <fr c i 
must be the vacuum solution for D > 2. In the case D = 2, we have E 2 = 0, so that 
((^cOi) 2 + ((^ci^) 2 = a 2 for all x. This requires that 3> c i (as |x| — * oo) has zero 
winding number and hence is in Q^°\ We can prove this result as follows. Let r, 9 
denote polar coordinates in the plane. For any nonzero r, <fr c i defines a map of the 
circle (with coordinate 6) to a circle (because of the condition (( < i ) ci)i) 2 + (( < I ) ci)2) 2 = 
a 2 ). The winding number n of this map cannot depend on r, as changing r is a 
continuous change. When r — > 0, all values of 9 represent the same spatial point, 
so that n — > 0. Hence n is identically zero which proves the result. 

2. Simplest nonlinear topological model 

The simplest model which modifies space topology can be found in one-dimensio- 
nal field theory. Consider the set of all mappings a from the real line R 1 (argument 
space) onto the circle S 1 (function space). S 1 can be parametrized by two real 
variables $ = ($i,$2), whose squares add up to one: $ 2 + <£ >2 l = 1- Note that func- 
tion space is not flat (circle S* 1 ). To prevent the escape of interesting structures at 
infinity we consider only the class of functions on S 1 , such that 3?(oo) = <&(— oo). 
The restriction of functions class allows us to identify argument space (line B}) 
with a circle. In other words it makes possible compactification of R 1 to a circleQ 

3 We h ave in mind confi gurations with (^"io) 2 + ('J'^) 2 = a 2 . Indeed, using Gauss the- 
orem from (1.1.6) and d[.l.^ ) we obtain / J* op d 2 x ~ / divPd 2 x ~ /^^(P • <K) ~ 0, if 



limixi^^ $ r (ii, X2, t) = 0. 

4 One point compactification theorem |5s| | (p. 86). Two-dimensional analog of the compacti- 
fication is known as a stereographic projection. 
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Figure 1.3. Example of one-dimensional model map in a simple 
nonlinear model. 



A field from n = 1 sector may be illustrated pictorially by a strip familiar with 
Mobius strip, except that the Mobius strip has a twist through tt whereas the sector 
n = 1 field has a twist through 2tt. ^^(x) then specifies twist angle of the strip 
about its center line at a given point x. Note quite different meaning we give to 
the external circle in this (see Fig. L3) and linear case (Fig. L2 y ). Since by classical 
field we mean a field that is single-valued under the action of the rotation group, 
it follows that particles involved must be bosons. Quantization of such a classical 
field introduces a quantity that can be interpreted as a particle number. It is known 
that after quantization the states, corresponding to n = 1 classical configuration, 
are, in fact, fermion states. Dynamics and quantum mechanical operators can be 
introduced^] into this theory [||, |59| |. 

To summarize, the simple nonlinear model fields are subject to nonlinear con- 
straint, when linear a model fields are not. This explains the need to reduce one 
argument space dimension in order to have the same global symmetry group for 
both models. In other words, function space of the first model is a flat space (vector 
space as well) , when function space of the second one is a compact manifold (not a 
vector space at all). 



3. Nonlinear a model in two spatial dimensions 

In this section we formulate nonlinear a model in fundamental SU(2) repre- 
sentation using well known Pauli and rotation matrices technique. The Skyrme 
model then is obtained by adding forth order term in field functions which ensures 
stable soliton solution in 3D. As a consequence, Skyrme model inherits all essential 
features from nonlinear a model. 

3.1. Formulation. It was sufficient to look at Q°° (the space of physical fields 
at spatial infinity) for the topological consideration in the linear a model. For 
solitons in nonlinear models it is often necessary to consider the topology of physical 
fields defined over all space. Physical fields in these models for all points x take 
values in a manifold J$t which generally is not a vector space. 



5 Lagrangian of the toy model: S£ = d^&sd^Qs', A 1 = 0, 1; s = 1, 2, where <l>, in addition, 
is subject to the constraint $| + $2 = 
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By definition group <&> acts on manifold Jtt transitively, if for any pair p, p' £ JH 
there exists an element g £ €>, such that T g {p) = p' . Assume, this is the case. Then 
Jttt is called a homogeneous space for If g p is the stability group of point p £ Jttt: 

g p = {he ®\hp = p}, (1.3.1) 

and Jtft is a homogeneous space for © then any two g p ,g p ' p,p' £ Jit are iso- 
morphic. If p 7^ p' and T g (p) = p' the isomorphism g p — > g p / can be defined^ 
h i — ^ ghg~ x . Now we can identify Jit with space of left cosets ©/jj by the following 
procedure. First let's fix point po £ JH. With each class of left cosets {<?g} we 
identify point T g (p ), where g = g Po is a stability group of point po £ Jtt. The 
identification is in one-to-one correspondence and do not depend on particular g in 
the class of left cosets. 

The nonlinear a model in two spatial dimensions has <&> =SU(2), g =U(1) 
= {e Mji , < a < 2ir} and J$t=6>/g is a two sphere S 2 . To show this, define 

K = gj z g f =n a ,J a ,, a' = (x, y, z) = (1, 2, 3), (1.3.2) 

where g £ SU(2) and J a > are generators of SU(2). 

K is an invariant under transformation^ 7r : g — > ge iaJz . A map 7r projects 
ge iaJ * £ for all a to the same point of left cosets {gg}. Since a is a contin- 
uous parameter, {gg} defines a two dimensional manifold, namely a two sphere 
SU(2)/U(1) = S 2 . Indeed, the scalar product of n a > is 1, 

TtK 2 = n a ,n a , = (n a ,) 2 = 1. (1.3.3) 

Fields $ a '(xi, x 2 ,t) of nonlinear 2-dimensional a model are subjects to the con- 
straint: $ a '( x , t)^a'( x , t) = 1 and thus can be identified with n a '. The action of 
on these fields yields 

K>=R*'v{9)*v, (I-3-4) 

where R a 'b' is the usual rotation matrix — an element of adjoint representation of 
SU(2): R £ 50(3). Note that the constraint is invariant under this action of (H>. 
The Lagrangian density is chosen so that it is invariant under 

Sf = -Zdp$ a r&>$ a ., [3 = const, = 0,1,2. (1.3.5) 



Because of the constraint on <& a ', the Lagrangian ( 1.3.5 ) does not describe a free 
system. Interactions of the field with itself are implicit. To see this, one can write 
in terms of two independent degrees of freedom, say 4>i = & x and $2 = & y 

( 2 , ) ( d ^k) 2 , $3^0, k is summed over 1 and 2. (1.3.6) 



2 Vl-$fc$fc, 

Note that the action of <& is nonlinear in terms of two degrees of freedom. 



6 It is assumed that group acts on from the left. Homogeneous spaces play an im- 
portant role in ensuring uniqueness of solution of equations of motion. There is no homogeneous 
space problem in 3-dimensional cr (Skyrme) model, because in this case there exists one-to-one 
correspondence between function space and SI.' (2) manifold, which, of course, is natural homoge- 
neous space for the same group. The problem again arises in SU(3) Skyrme model, when SU(2) 
ansatz is employed. 

7 And only under these transformations, because [J z , J x [y]] 7^ 0- 
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Figure 1.4. Example of mappings: S 2 sphere — > S 2 sphere. 



The energy density associated with ( I.3.5| ) is 



<?(*,t) = ^((d <i> a ,) 2 + (d^ a ,) 2 ); i = i,2. 



(1.3.7) 



3.2. Topological structure. The vacuum solution (which is subject to the 
constraint on <fr) is 3> V ac = const. S is invariant under global SU(2) transformations. 
Thus 3>vac can be reduced to (0,0,1) by action of SU(2) without affecting the 
energy {$ = for <fr V ac)- After the choice only rotations about the third axis 
leave the vacuum invariant, consequently, the global SU(2) is spontaneously broken 
to a global U(l), by special vacuum choice. Next, consider general configuration 
with nonzero, but finite energy. Let us first show that homotopic sector generation 
mechanism described in linear a model in Sec. [l] now fails. Indeed, for large |x| = R 
the fields 3? define a mapping from circle S 1 with radius R to S 2 . Since 7Ti(5' 2 ) = 
this mapping is homotopic^ to 3? V ac at |x| = R. Despite the fail there is a way out. 
Indeed, finiteness of energy requires 3>(x } i) to approach <fr vac as r — > oo and that 
the rate of approach is fast enough to guarantee that the energy E is finite. Again 
we choose 3? V ac = (0,0,1). Assume that 4> approaches 4> V ac and there is no angle 
dependent limit at r = oo. Thus, we may think of all points at spatial infinity as 
being a single point. Such a restriction of function class 3> essentially converts (in 
topological but not metrical sense) the plane R 2 = {(ari,ara)} at a constant time to 
to the surface of a two sphere S 2 . The fields $ = ($i, $ 2 , $3) are well defined on 
this S 2 in view of boundary conditions. Consider this in detail. Let £ a , a = 1, 2, 3 



8 Even if the fundamental group of the manifold is nontrivial: 7ri(jH) ^ 0, the fields at 
R would have to be homotopic to the vacuum solution $ vac . This is so since # defines a trivial 
mapping at r = 0, and the topological index cannot change as r is continuously varied from r = 
to R. (Constraint 3> 2 = 1 here is fulfilled at each point. In linear a model, in contrast, we had 
this condition satisfied at infinity only.) 
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Figure 1.5. Examples of mappings: cylinder — > cylinder and 
cylinder — > Mobius strip. 



be the stereographic coordinates associated with x 

^( X ) = ^Y> « =1,2; 6(x)-^j 

£ 2 (x) + £ 2 (x) + £ 2 (x) = 1, r 2 = x Q x Q . (1.3.8) 

Coordinates £ span a two-sphere S 2 . They are not global valid coordinates for R 2 
which unlike S 2 is not a compact manifold. Indeed, all the "points of infinity" of 
R 2 which correspond to r — > oo are mapped to one point of S 2 , namely (0,0,1). 
In reality the "points of infinity" are not points at R 2 at all. Thus, to get a 
topologically accurate representation of R 2 we should remove the north pole (0, 0, 1) 
fromS 2 : R 2 = S 2 \{{0, 0, 1)}. 

The topological difference between R 2 and S 2 can make difference for some 
functions. For example, the function /(x) = |x| is a continuous function on R 2 , 
but the function obtained by the substitution x a = £ a /(l — £3) is not a continuous 
function on S 2 , becoming infinite at north pole £3 = 1. Another example is the 
function x = x/|x| which is continuous on R 2 , while its image function on S 2 has 
no well defined limit as the north pole is approached. However, for functions which 
approach a constant limit as r — > 00, the change of variable x — > £ does produce a 
well defined function on S 2 . In this sense then, because of the boundary condition 
on we can imagine that the space on which the field <fr is defined is S 2 . 

Thus, the configuration space Q of the nonlinear a model is made up of fields 



which map S 2 to Jtt = S 2 * : S 2 -> S 2 (see Fig. |D]). The situation, thus, is 
analogous (despite quite different origin of the map) to the maps $00 in linear a 
model, where we had S 1 — ► S . It is then plausible to expect that for nonlinear 
a model the configuration space Q falls into an infinite number of disconnected 
components Q n , with Q — U„ Q n . This result is true. Here n is generalization of 
the previous winding number associated with <froo and is also called the winding 
number. It indicates the number of times the sphere S 2 is covered by sphere S 2 as 
£ runs over all values. Strictly speaking, one cannot transfer this very illustrative 
definition from the one-dimensional case (number of times one circle is covered by 
another circle) to higher dimensions (number of times one S k sphere is covered by 
another S k sphere) . Mathematically irreproachable definition of mapping degree^] n 

9 The illustrative definition of the winding number is indeed correct for the one-dimensional 
case. The reason is that curve does not have internal structure. In general it is not true for higher 
dimensional manifolds. For example mappings cylinder — > cylinder and Mobius strip — » cylinder 
in Fig. |l."i| obviously are not homotopic, but illustrative winding number definition does not allow 
us to clearly distinguish both cases. 
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on m > 2 dimensional manifold is based on very general triangulation concept |60|. 
Despite the criticism, it is known that for mappings S — ► S k , k = 1, 2, . . . the k- 
th homotopy group is nk(S ) — Z, which in some sense justifies the illustrative idea 
of the winding number. After the generalization once again, equivalence classes Q n 
can be made into a group under a suitable product. This group is called the second 
homotopy group and is denoted by 7r 2 (Jttt). Here JH = S 2 . Like 7Ti(S' 1 ), tts(S ) 
is isomorphic to the group Z of the all integers under addition. 

The equivalence class Qo contains the vacuum solution * (0) = * vac = (0, 0, 1). 
Qo consists of all maps which are homotopic to 3> V ac- An element of Q\ is ob- 
tained by simply setting t) = here t being fixed and £ being coordinates 



defined in flL3.8|). A typical element of Q [n) is 



cosnip, £i (x) = sin -d cos ip, < $ < tt, (1. 3. 9a) 

<J> 2 n) (x,i) = sintfsinn^, £ 2 (x) = sin •& sin </>, < ip < 2ir (I.3.9b) 
$3 n) (x,t) = cosz? £ 3 (x)=costf . (I.3.9c) 

Here y> are spherical coordinates for the argument two-sphere S 2 . Q n consist of 
all maps homotopic to the & n \ 

Again, the significance of the above classification (since time evolution is a 
continuous operation), is that the integer n is an integral of motion. It is useful 
to have an explicit formula for this conserved quantum number. For this purpose 
consider the current 

j"(x,t) = -^w e! a ' x * a ,d v *b,dx*c, (1.3.10) 

where e^ uX is the totally antisymmetric tensor. This current is conserved regardless 
of the equations of motion. Taking its divergence we obtain 

d^i^t) = ~e a , v& e^%<S> a , d v $ v d x <S> d . (1.3.11) 

07T 



The right hand side of ( 1.3.11 ) contains the triple scalar product of the three tan- 
gent vectors do*, and 9 2 3? defined at (x, t). When multiplied by d 2 x.dt, it 
represents an infinitesimal volume element at (x, t) . But because of the constraint 
on <&, the tangent vectors at (x, t) are enforced to lie in a plane. Conseq uently , the 

we 



dj° 

have — — = 0. It follows that the associated charge^ 



volume element and the right hand side vanishes <9 M j M = and from ([L3.11 
df 
dt 



= / d 2 xe rt . c - ey&a'd&vdj&cr, (1.3.12) 



is a constant of motion. Its value is the conserved quantum number; it has the 
value n when <fr = S Q n . The factor — ^ is chosen so that B(&) is, i n fact , 



an integer. To see that B(&) = n, write $ in terms of spherical coordinates ( [.3.9 ). 
Then 

= ^- J frini?(x)dp(x)di?(x). (1.3.13) 



10 Note that the charge density now is not a pure divergence j° 7^ divP (as it was in linear 
model) and has only integer values when properly normalized. 
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Since ■^sm.'&dip A di) is the normalized volume element^ on two sphere, B($>) 
indicates the number of times the sphere S 2 is covered as x runs over all values 
and is, therefore, an integer. Derrick scaling argument rules out (see Sec. 1.3 ) 
the possibility of having nontrivial static solution to a linear scalar field in two (or 
greate r) space dimensions. However, for the nonlinear a model with the Lagrangian 
( t.3.5 ), Derrick's argument can only be used to rule out the existence of static 
solutions in all but two space dimensions. This is because the static energy contains 
only one term which we denote by E a . Under x — > Ax, it scales like E a — > X 2 ~ D E cr . 
The minimum value of the energy for this variation is zero in all except D = 2 
dimensions. 

A lower bound on the energy (the "Bogomol'nyi bound") [54, [55| for the 
classical solutions can be obtained from the identity 

(d i $ a ,±e a , b , c ,e ij $ b ,d j $ c ,f >0. 

After completing the square, we can write 



(1.3.14) 



~E S = J d 2 x($*«') 2 > 87r|n|. 



Here the (3 is the same as in ( [.3.5 ). The bound is saturated if 
Here identities 



(1.3.15) 



(1.3.16) 



Sfc" <W' and Cij e< 



(1.3.17) 



Ca'b'c' £a'b"c" — h'b" 5c< 

have been used. 

A general solution to equation ( 1.3.16 ) was obtained by A. A. Belavin and 
A.M. Polyakov [61|. Here we shall only look for a spherically symmetric n = 1 
solution. Spherical symmetry in two spatial dimensions means eijXidj$ a ' = 0. 
This condition is consistent with the constraint on However, it has the unde- 
sired result that all fields satisfying it have B(&) = 0. This is because the general 
solution to eijXidj& a ' = is 



$ a '(x,i) = $a'(r,*), so that = x. 



(1.3.18) 



or r 

Upon substituting this into (|l.3.12| ) we immediately obtain the result B(G>) = 0. Af- 
ter some symmetry requirement modification it is possible to obtain configurations 
with ^ 0. We refer for details to @ (p.119-121). 



3.3. Going to 3D space. The change of space dimension is a highly non- 
trivial action. The existence of many objects and phenomena which are allowed in 
some dimensions are forbidden in another's. For example, two-dimensional crea- 
tures should have different digestive tract and blood circulation system, otherwise 
eating or blood circulation would divide them in two separate halves 1 62 1 (p. 164). 
There also would be problems with more than three space dimensions, in particular 
with gravitational force. As a consequence, orbits of planets would be unstable. 
Here are theories more or less successfully describing phenomena when higher di- 
mensions are introduced (string theories). The problem usually then becomes how 
to reduce these nonobservable dimensions. Our aim now is to construct realistic 3D 
nonlinear theory, with essential features inherited from two-dimensional a model. 



llr The symbol A is an exterior multiplication mark. 
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What do we need in order to extend a model to real 3D space? First, we note 
that the compactification method used in nonlinear a model can be easily extended 
to the 3D case. Indeed, compactification R 3 at fixed time t to S 3 leads to the 
mapping S 3 — > S 2 with trivial homotopy group 7T3(5 2 ) = 0. Therefore in order to 
get nontrivial topological classes we should add one more field, satisfying 

*;(x,t) + $|(x,t)+*|(x,i)+*J(x,t) = l; x= (x 1 ,x 2 ,x 3 ). (1.3.19) 

Then again we have 7r 3 (5 3 ) = TL. Additional field component ensures that field 
<fr has values in the whole SU(2) manifold. Thus group manifold becomes natural 
homogeneous space for the group itself and no identification of ill with space of 
cosets is needed. The additional field, however, does not eliminate the soliton 
stability problem. The simplest way to eliminate Derrick scaling argument (which 
excludes static stable nontrivial solution) in classical level of the theory^ is to add 



LU 



a new term in the Lagrangian density which would stabilize the energy (I 

Skyrme succeeded in suggesting the following four th-or der term (ensuring stable 
soliton solution) to be added to Lagrangian density ( I.3.5| ): 



&Sk =^Tr[R /1 ,R I/ ] [R",R"] 

^TrfRo.R] [R°,R] +^Tr[R-,R J ] [R\R J ], (1.3.20) 



where 



R M = (S p U)Ut; U(x,i) e 5*7(2). 

The contribution to static soliton energy comming from the Skyrme term scales as 

Esk(l) - E Sk (X) = X 4 -° E Sk , (1.3.21) 

under a scaling transformation U(x) — * U(Ax). Requiring again that A = 1 corre- 
sponds to energy minimum yields the equation 

E a = E sk , (1.3.22) 

in D = 3 space dimensions. Assuming that soliton energy is proportional to its size 
R and taking into account dimensions of /^[MeV] and e[dimensionless] we conclude 
that the leading term (E a ) is proportional to ~ af%R, whereas the Skyrme term 
t o ~ (ci,C2 are positive constants). Thus, adjusting soliton size R equation 
( t.3.22 ) can always be satisfied^, for nonzero soliton size R. 

We can also add terms involving more than four derivatives (for example, 
and Jzf 6 ' terms in Sec. ^|.|l.5|). There is no good argument to suggest that these 



terms are ignorable. For example, the so-called large N c limit of QCD |6S , 70] fails 
to show that higher derivative terms are down by powers of N c as compared to the 
leading terms. Despite these criticisms, we will approximate the action density by 

J?? CT + Sfsk ■ 

When the Skyrme term is also included in the Lagrangian density, there is 
an elegant lower bound to the energy of soliton. The bound is analogue of the 



12 Ther a exist stable solutions with only a term when coupling to vector mesons is in- 
cluded fl63| . There is discussion on the market, however, whether scale parameter or breath- 
ing mode quan tization can stabilize the solution wi tho u t th e Skyrme term. For arguments 
see |p4| , |65| , |66fl , for con tra-arguments we refer to |3l], |67| , |68|. 

^^Vloreover, (f.3.22| ) is satisfied for only one positive -R value due to the second-order algebraic 
equation ci/^e 2 i? 2 — C2 = 0, describing the energy extremum condition. 
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Bogomol'nyi bound we considered earlier although it predates the latter by many 
years. It is based on the observation that 

J d 3 xTr(^R, + le^fR^Rfe]) 2 < 0. (1.3.23) 

To show this result, notice that and e^-fc [R 3 - , R*] are antihermitian matrices 
and that for any antihermitean matrix G TrG 2 < 0. From this we arrive at the 
bound 



1 



,]"xTrl — R 

1 4 ' ' 32e2 



' ° 2 ' [Ri,Ri] 2 ) > ^^\B(V)\, 



(1.3.24) 



due to Skyrme. Here |S(U)| denotes a wind ing number, explicit exp ression of 
which is given in the next chapter (see ( [1.1.18 )). The left hand side of ( L3.24| ) is 
the potential energy of field U(x). The bound thus shows that in the presence of 
the Skyrme term, the soliton energy and mass are bounde d from below. Although 
there is no nontrivial solution which saturates the bound ( 1.3.24 ) , static solutions 
to field equations are known to exist for \B\ = 1. In Sec. [2 of Chapter || we shall 
discuss the "spherically symmetric" static B = 1 solution, which by ~ 23% 
exceeds the boundPL 



4. QCD and the Skyrme model 

Links between fundamental theory (QCD) and phenomenological theories of 
strong interactions (including the Skyrme model) are briefly considered here. 



4.1. Historical remarks. T.H.R. Skyrme proposed his model in 1961 |71|. 

For almost two decades the theory has been ignored and only in early 80-ies it 
has been realised that the model, as effective theory of mesons, may provide a link 
between QCD and the familiar picture of baryons interacting via meson exchange. 
Low energy domain of QCD becomes forbiddingly difficult due to the rising cou- 
pling constant which possess a major obstacle to a satisfactory description of the 
dynamical behaviour of the elementary quark and gluon fields of QCD at the rele- 



vant large distances. R. Rajaraman's [54| and E. Witterfs |69| results suggest that 



baryons may be regarded as soliton solutions of the effective meson theory without 
any reference to their quark content. This was precisely what Skyrme had sug- 
gested in his remarkable papers JH], fn], 72 1. There are a lot of works analyzing 



one or another aspect of this extremely important and interesting problem. For 
overview we refer to Ref. |flo|| and references therein. Here we consider only gen- 
eral phenomenological requirements for effective theory of strong interactions and 
very briefly describe the 1/N C expansion idea. Unfortunately, we completely escape 
chiral perturbation theory recently making a huge progress. This theory, however, 
explicitly involves baryon fields (when describing processes involving baryons) and 
is outside the Skyrme's idea that baryons are solitons of meson fields. 



14 We shall see in Chapter ^ that the negative quantum mass correction (HI. 2. 55) can lower 
quantum Skvrmio n mas s. The question, however, can be asked, whether quantum Bogomol'nyi 
bound similar to (1.3.24) can be defined when dynamical variables don't commute. 



4. QCD AND THE SKYRME MODEL 



13 



4.2. General requirements for effective theory of strong interactions. 

The starting point is an idealized world where Nf\ avouIS = 2 or 3 of the quarks are 
massless (u, d and possibly s). In chiral limit the QCD Lagrangian exhibits a global 
symmetry 

SU(JV»£ x SU(JV»fl xU(l)y x U(1) A . (1-4.1) 

chiral group G 

At the effective hadronic level the quark number symmetry U(l)y is realized as 
baryon number. The axial U(l)^ is not a symmetry at the quantum level due to 



the Abelian anomaly |73, 74 1 that leads, for instance, to M v ' ^ even in the chiral 
limit. 

There is compelling evidence both from phenomenology and from theory that 
the chiral group G is spontaneously broken [ f75| |: 

• Absence of parity doublets in the hadron spectrum. 

• The Nj — 1 pseudoscalar mesons are by far the lightest hadrons. 

• The vector and axial- vector spectral functions are quite different. 

• In vector-like gauge theories like QCD (with the vacuum angle #qcd)> vector 
symmetries like the diagonal subgroup of G, SU(7V/)y, remain unbroken. 

All these arguments together suggest very strongly that the chiral symmetry G is 
spontaneously broken to the vector subgroup SU(iV/)y (isospin for Nf = 2, flavour 
SU(3) for N f =3) 

G -> H = SU(N f ) v . (1.4.2) 

Then the Goldstone theorem tells us that there exist Nf — 1 massless mesons. 
For two flavours, these Goldstone modes are identified with the three pions, while 
for three flavours, these modes are identified with the pseudoscalar octet. In chiral 
limit (when quarks have zero masses), the pseudoscalar mesons are exactly massless. 
They become massive when the interactions between the quark and Higgs fields are 
turned on, the quarks acquire mass and G gets explicitly broken in the Lagrangian. 

The effective Lagrangian emerges when we attempt to construct a model which 
describes the dynamics of these Goldstone modes. Let us list the properties we 



require for this Lagrangian in the zero quark mass limit |51|: 

1. The Lagrangian L must be invariant under G = SU(JV/)l x S\J(Nf)n, this 
property being the analogue of the G-invariance of the QCD Lagrangian. 
Thus L is to be constructed from a multicomponent field 3> which is trans- 
formed by G, L being invariant under these transformations. 

2. Field should have exactly Nj — 1 degrees of freedom per space-time point. 
This is a requirement of minimality: we want to describe the dynamics of 
the Goldstone modes and only of these modes. It is possible to improve 



effective theory by introducing vector or /and axial vector mesons [24|, or 



even massive non- Abelian gauge bosons | |76| |. 
3. We require that the subgroup of G which leaves any value of the field invari- 
ant is exactly (or isomorphic to) subgroup H and no more. If this can be 
arranged, then we would have nicely built in spontaneous symmetry break- 
down G — ► H in the geometry of the fields itself. 

It is an easy task to check that the Skyrme model satisfies all these require- 



ments J44j 



14 



I. INTRODUCTION 




k 



i 



g/vm g/vm 



Figure 1.6. Gluon correction to quark propagator in standard (a) 
and 't Hooft-Witten (b) notations. 

Chiral perturbation theory (CHPT) is also based on similar chiral symmetry 
principles. Generally, many more chiral invariant terms can be included into the 
Lagrangian. The Skyrme model, basically, takes only two of them L 2 and L4. 
CHPT, on the other hand, provides us with a scheme which tells us which terms 
should be included and which ones should not. Roughly speaking, the essential idea 
of chiral perturbation theory is to realize that at low energies the dynamics should 
be controlled by the lightest particles, the pions, and the symmetries of QCD. 
Therefore, S-matrix elements, i.e. scattering amplitudes, should be expandable 
in Taylor-series of the pion momenta and masses^, which is also consistent with 
chiral symmetry. This scheme is valid until one encounters a resonance, such as 
the p-meson, which corresponds to a singularity of the S-matrix. It should be 
stressed, however, that chiral perturbation theory is not a perturbation theory in 
the usual sense, i.e. it is not a perturbation theory in the QCD coupling constant. 
In this respect, it is actually a nonpertubative method, since it takes infinitely 
many orders of the QCD coupling constant in order to generate a pion. In the 
meson sector CHPT is quite successful, whereas the precision achieved in heavy 
baryon CHPT is not comparable to the meson sector accuracy. For explanation we 
refer to lectures ||75|. 



4.3. The 1/N C expansion. Assuming confinement, the asymptotic states of 
QCD are not the coloured quarks and gluons, but rather the observed colour singlet 
hadrons. In view of this, one might wonder whether in some way QCD itself could 
not be equivalently formulated in terms of these observed asymptotic degrees of 



freedom. Quite remarkably, the work of G. 't Hooft Hrt)H and E. Witten [|69| shows 
that QCD is indeed equivalent — in the full field theory sense — to a theory of 
mesons and glueballs^, with meson-meson coupling constant ~ From the 

first sight, there seems to be one very large gap in the equivalence 

confined QCD = theory of mesons and glueballs, 

namely, where are the baryons? It is here that the real interest of the 1/N C idea 
lies. E. Witten showed |69[ that for large N c baryon masses scale like N c . This is 
reminiscent of the behaviour of solitons in a theory in which the coupling constant 
is g: the soliton mass is ~ 1/ff 2 , so that putting g ~ we find mass ~ N c . 

But this interpretation is exactly what Skyrme suggested. 



15 In the baryon sector one has an additional parameter — the nucleon mass. 
16 There exist at least few effective field theories in four dimensions, as the number N of fields 
of some type becomes large |77||. 



4. QCD AND THE SKYRME MODEL 



15 




Figure 1.7. Gluon loop in standard (a) and 't Hooft-Witten (b) notations. 

Let us now briefly explain why meson-meson coupling constant should scale like 
~ l/N c . To this end let us reformulate QCD for an arbitrary number of colours N c . 
In such a theory there are N c quark degrees of freedom and N% — 1 « AT? gluonic 
degrees of freedom (for large N c ), Consider then the simple gluonic correction to 



the quark propagator depicted in Fig. 1.6 a 



Even after we specify the colour index of the external quark, this diagram 
receives a combinatorial factor of iV c corresponding to the N c possible values for 
the index of the internal quark. This is easy to see in 't Hooft-Witten notations 



Fig. 1.6,, where gluon in combinatorial sense (and only in this sense) is equivalent 
to quark-antiquark pair. The resulting loop corresponds to the summation over all 
possible quark index values and is responsible for the combinatorial factor N c . 

If we want the theory to have a smooth — but nontrivial — limit as N c — > oo we 
must compensate this combinatorial factor. Thus we require that the qqg vertex 
scale like ~ l/y/N c . The same result is obtained for the trilinear meson-meson 
(gluon-gluon) coupling constant as can be seen from diagrams Fig. [L7| a and Fig. 1.7 



G. 't Hooft noted that not all diagrams are of the same significance when 
N c — > oo. Simple power counting, similar to those just described, implies ||70) 
that in this limit only planar diagrams become important. Analysis of all planar 



diagrams [ 70 ] , which is out of the scope of this work, together with confinement 
assumption of large N c QCD leads to the following conclusions: 

• Mesons are stable and to leading order non-interacting particles. Their 
number is finite. 

• Amplitudes of elastic meson-meson scattering are of the order of ~ l/N c 
and are expressed as sum only of tree level diagrams]^]. 

• Baryon- meson scattering can also be analysed in similar fashion [12, |l3|, 



78|, |7g|, |80|. 



In other words, QCD seems to reduce smoothly to an effective theory of mesons 
(and glueballs) with the effective coupling constant of the order of ~ l/N c . 



Tree level diagrams in this case describe only meson (but not gluon or quark) exchange. 



INTRODUCTION 



CHAPTER II 



Classical Skyrme model 



The chapter deals with classical Skyrme model. Classical model assumes that 
dynamical model variables and its time derivatives commute. The model is usually 
formulated in rotation Rij and Pauli r matrix formalism, illustrated in 2D nonlinear 
a model. Wigner matrices and su(2) algebra generators J represented in circular 
basis are, nevertheless, more convenient for model formulation in arbitrary reducible 
representation and, therefore, will be followed in this and subsequent chapters. 

1. Formulation 

The section serves as a formulation of the classical SU(2) Skyrme model in 
arbitrary irreducible representation. An emphasis is put on expression dependence 
on representation. Physical quantities (mass, coupling constants, etc.) are inde- 
pendent of the representation after the proper model parameter renormalization is 
employed. 

1.1. Parametrization of the symmetry group. Chiral group is a group of 
transformations in the internal (isotopic) space, under action of which left and right 
states transform independently. Simple nonabelian (six-parameter) chiral group 
is obtained by multiplying two rotation groups directly SU{2)i, B , i ® SU{2)n ight . 
There is no linear realization of the group in 3D isospace. One can choose either to 
extend the isospace to 4D, where the linear representation exists or to construct the 
nonlinear representation. The natural nonlinear representation (which we follow 
further in the work) is obtained when the group parameters space (manifold) is 
identified with the space where the abstract group transformations are realized. 
For example, SU (2) matrix in well-known Euler-Rodrigues parametrization ao,a 
takes a formQ j8l| 

U(oo,a)= ( ao + m3 mi+a2 )=a -l + l a-T, (II.l.l) 

where the group parameters spaceQ ao,a itself is restricted by the constraint 

o + a 2 = l; a ,aelR. (II.1.2) 



a, 



The presence of constraint ( II, 1.2] ) gives rise to additional problems in quant: 



zation of the theory. From our point of view unconstrained parameters are more 



suitable for this purpose, but see |43| . Such an unconstrained parameters are, for 
example, triple of Euler angles a = (a 1 , a 2 , a 3 ) |82[ 



< a 1 < 2tt, < a 2 < tt, < a 3 < 4tt. (II.1.3) 



1 T.H.R. Skyrme |7l| formulated his model in terms of a and pion 7r-fields U = (cr + ix- 7r). 

2 Sometim es the set ag,a is called a 4-isovector. The reader should be aware that because of 



the constraint (II. 1.2) the set isn't a vector space. 
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An arbitrary reducible SU(2) matrix in the Euler angles parametrization can be 
expressed as a direct sum of Wigner D J (o;), (j = |, 1, |, . . . ) functions. 

1.2. The Lagrangian. For reasons of simplicity and without lose of general- 
ity0 let us formulate the model in the arbitrary irreducible SU(2) representation j. 
Euler angles c*(r,t) become the functions of space-time point (r,i) and form the 
dynamical variables of the theory. Model is formulated in terms of unitary field 

U(r,t) = D J (a(r,i)), (II.1.4) 

all physical quantities being functions of this field U. In the quark picture the 
analogue of is the complex 2x2 matrix q 1 1 ~ 75 q\ corresponding to pseudoscalar 
mesons ||53|| . Note that this analogue is only valid in the fundamental representation 
of SU(2) . Unitary field can also be expressed in terms of pion fields 7r and unphysical 
<x field 

U(x) = — (cr + ir-n). (II.1.5) 

The basic Skyrme model is described by chirally symmetric Lagrangian density^ 

f 2 1 
^ Tr{R M R^} + — 

4 ' 32e 

where the "right" current^ R^, known for mathematicians as Maurer-Cartan form, 
is defined as 

R M = (9 M U)Ut, (II.1.7) 

f-„ (pion dec ay con stant) and e being parameters^ of the theory. Let us explore the 
Lagrangian ( Il.l.tj ) algebraic structure more closely. To this end it is convenient to 



S£ = -if Tr{R M R M } + ^ Tr{[R M , R„] [R M , R"]}, (II.1.6) 



introduce a contravariant circular coordinate system. The unit vector x in these 
(contravariant) circular coordinates is defined in respect to Cartesian, spherical and 
circular covariant coordinate systems as 

x +l =- —=(xi - ix 2 ) =-^sini?e- lv = - X- X , (II.1.8a) 

x° ~X3 = cos$ =xo, (II. 1.8b) 

x~ 1 = —={x\ + ix-i) =—j=sni{)^' ip = — jc+i, (II. 1.8c) 

v2 v 2 

respectively. Then the general inner (scalar) product of two algebra elements can 
be defined as 

Ti(j-\J a J b \j-) = {-l) a \j(j + l)(2j + l)S a .- b , (11.1.9) 



3 



Formulation of the model in arbitrary redi 



cible representation simply involves summation 
over representation j and is explained in Sec. p. 4[ 

4 We consider chiral transformations in detail in Sec. Q and Sec. y of Chapter 

5 The theory as well can be formulated in terms of "left" current L M = U'(8 M U). Recent 
political tendency rendered "right" current more popular, what we believe defined our choice. 

6 Note, however, that the parameter f n value cannot be determined within the framework 
of strong interactions only, because pions are by far the lightest strongly interacting particles 
and, thus, are stable in this theory. Experimental value of f % is 93 MeV. It is claimed | |7o| | 
that par ameter e value can be extracted from the tttt scattering data using formulas given in 
Ref. | p3| . The result is e = 7.4. The Skyrme constant e also has been roughly estimated by 
assuming that the Skyrme term arises by "integrating out" the effects of a p(770) meson; this 
yields e = m p /(2f^) = 5.83 E3. 
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where the su(2) generators J satisfy the commutation relation 
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[Ja, Jb 



1 1 1 

a b c. 



Jc 



c = a + b. 



(II.1.10) 



The factor on the r.h.s. in ([1. 1.10) is the Clebsch-Gordan coefficient (la 16|lc) in 
a more transparent notation. Wigner D J function parametrization in the form 

(j-\W(a)\j-) = (j-|exp(zV2Vj ) exp(-a 2 ( J+ + J_)) cxp(iV2 a 3 J ) \j-), 

(IL1.11) 

makes it easy to obtain the following relations: 



A Dl n {oc) = C\ a \ a ) (jm\J a \jm') D> m , n (a), 
A D* mn (-a) = -C[ a \ a )Dl nl (~ a ) (jn' \J a \jn) , 
A Dl n {cc) = C'^( a )Di nm ,( a ) (jm' \J a \jn) , 



t) 



^- Di m (- a ) = -C^ a \cx) (jm \Ja\jn) l£, n (-a), 



where the coefficients 

Cl a \a.)=Dl al {<x)C'} a ' >(a) 
^ (o) (a)=£»i, ,(-a)C^ ' ) (a) 



have the explicit form [49| 



C| +) (a)=0 C^(a)=-e- M " C^a^ -i sin a* e 



.(+) 



a = (ai,a 2 ,Q!3), 
-a ee (-a 3 , -a 2 , — ai), 

(+)/ 



(t\a)=0 C ( -\ot)=- 



C { °\a) =iV2cosa 2 
C { s ~\a)=i sin a 2 e* Q 



and satisfy orthogonality relations 

£ C7 (m) (a)^ (a)= £ C? m > («)C* , (a)=<^ , 

m m 

^Ci m) (a)q n) (a) =^C7; (ro) (a)C ( ^(a) =5 m ,„. 



(11.1. 12a) 
(II.1.12b) 
(II.1.12c) 
(II.1.12d) 



(11,1. 13a) 
(II.1.13b) 



(II.1.14a) 
(ILL 14b) 
(II.l,14c) 



(11.1.15a) 
(II.1.15b) 



Using formulas (|LL12j) the right current R M can be reduced to the form 

(«), 



{Rn)mm' = d^Cf'ia) (jm\J a \jrri), 



(II.1.16) 



and clear ly have values in su(2) algebra. Relations (|LL9|), ( [I.L10| ) together with 
formula ( [I.1.16| ) allow us to express the Lagrangian density ( II.1.6 ) in terms of the 
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Euler angles [49 1 

2 =~j(j + l)(2j + 1) f^(d M aW + 2 cos a 2 d^a 1 O^a 3 ) 

( d„a 2 d»c?(d v a l d v a l + d v a 3 d u a 3 ) - (d^d^a 2 ) 2 
16e z V 



{d^d^a 3 ) 2 + sinV^a 1 ^ 1 d v e?d v u 3 - (d^d^a 3 ) 2 ) 



+ 2cosa 2 (d M a 2 <9 /i a 2 d v a x d v a 3 - d^d^d^d^a 3 



(III. 17) 



The only dependence on the dimension of the r eprese ntation is in the overall factor 
j{j + + 1) as it could be expected from ( II. 1 .Sj ) . This implies that the equa- 
tion of motion for the dynamical variable a is independent of the dimension of the 
representation. 

Note. We introduce additional normalization factor 1/N in the definition of quan- 
tum Skyrme Lagrangian in Chapter HI. The motivation comes from considerations 
below. 



1.3. The topological current. The fol lowin g con structi on of "right" cur- 
rents is called topological current density (cf. ( 1.1.5 ) and ( I.3.1C )): 

^^TrlR^RJ. 



3 • 8iV7r 2 



(II.1.18) 



The integral associated with (EL1.1S) is a conserved quantity. The normalization 
factor N depends on the dimension of the representation and has the value 1 in 
the fundamental (j = |) representation. The baryon number^ B is obtained as the 
spatial integral of the time component S3® . In terms of Euler angles a the baryon 
current density takes the form 



1 



3 • 6 • 87Vtt 2 



j(j + l)(2j + l)sina 2 e' il//37 e aa d v a i dp(J'd 1 o}. 0X1.19) 



As the dimensionality of the representation appears in this expression in the same 
overall factor as in the Lagrangian density ( II. 1.17 ) it follows that all calculated 
dynamical observables will be independent of th e dime nsion o f the re prese ntation 
at the classical level. The same overall factor in ( II.1.19| ) as in ( [1.1.17 ) and ( Il.l.gj ) 
also indicates that the topological (or baryon) current density 38^ can be expressed 
in terms of scalar product of algebra elements. This is indeed the case | 56 | 



e^Tr^R^RjRJ. 



(II.1.20) 



The forms (II. 1.18), ( II.1.2C| ) make no dif ference for ba ryon nu mber B in the quan- 
tum case. This can be seen both from ( II. 1.18 ) and ( [1.1.20 ) a s time derivatives 
Ro are not involved in the expressions. A more symmetric form ( [1.1.18 ) is usually 
used. 



7 Topological index (due to its conservation) is identified with baryon number in the Skyrme 
model. The following expressions are used as synonyms in the physical literature: topological 
index, Chern-Pontryagin index, winding number, soliton number, particle number, baryon number. 
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Figure II. 1. Hedgehog configuration. Arrows indicate the direc- 
tions of the isovector field F(r) at different points in coordinate 
space. 



1.4. The hedgehog ansatz. The general s olution of equations of motion, 
which follows from variation of the Lagrangian ( [1.1.17 ), is not found. Skyrme 
suggested the static soliton solution in the fundamental representation of SU(2) 



U 



= e i(T-f)f(r) _ 



(II.1.21) 

Here r is isovector of Pauli- isospin matrices and f denotes unit spatial vector Th e 
object described by ( [1.1.21 ) has a very peculiar geometric structure (see Fig. [LI): 
at each point x in 3D space the associated isovector F(r) points in a radial di- 
rection with respect to the spatial origin x = 0, where the centre of the object is 
located. This rad ial stru cture has prompted the handy name of "hedgehog" for 
the configuration ( [1.1.21 ). In order to find its generalizations for representations 

1 /2 

of higher dimension one may compare it to the matrix elements Dmn(ot), and thus 
obtain the explicit expressions for the Euler angles a in terms of the chiral angle 
F(r). The resultQis [|§ 

a 1 = ip - arctan(cost? tanF(r)) - tt/2, (H.1.22a) 

a 2 = -2 arcsin(sintf sinF(r)), (H.1.22b) 

a 3 = -<p- arctan(cosi? tanF(r)) + tt/2. (II. 1.22c) 

Here the angles <p, $ are the polar angles that define the direction of the unit vector 
f in spherical coordinates. 



Substitution of the expressions flll.1.22; ) into the general expression (jllTJ) for 
the unitary field U then gives the hedgehog field in a representation with arbitrary 



8 When yj,d, F(r) in ( p.l.22| ) run over values < ip <2ir, < < tt, 0< F(r) < it, the 
range of ai,a2,«3 is —it < a\ < 2tt, —n < cti < 0, — 2tt < 03 < tt and thus differs from 
([1. 1.3) range. This, however, can be fixed by dividing the parameters area in a proper way and 
moving each part by some fraction of tt. 
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j. As an example, the hedgehog field in the representation j = \ has the form 



G isinFsin^e"^ 
ismFsm-de^ G* 



(IL1.23) 



where we have used abbreviation G — cosF + isinFcos$. For j — I the same 
substitution yields 



G 



iV2smFsm$e- lv G -(smFsintfe"^) 2 



U 



iv^sinFsintfe^G 1 — 2 sin 2 F sin 2 $ iv^shi-Fsintfe^G* 
-(sinFsintfe^) 2 iv^sinFsintfe^G* G* 2 

(II.1.24) 



The Lagrangian density (II. 1.17) reduces to the following simple form, when 
the hedgehog ansatz ( II. 1.22 ) is employed: 



^ = --j(j+l)(2j+l) 



■ sin z F 



1 sin F , 



2F'- 



sin 2 F 



(11,1.25) 



For i = 5 this reduces to the result of Ref. 0], The corresponding mass density is 
obtained by reverting the sign of Jz? , as the hedgehog ansatz is a static solution. 

The requirement that the soliton mass be stationary yields the following equa- 
tion for the chiral angle F(r) pi]: 



2 , sin 2F 

~ F 



1/1 



■sin 2 F sin2F 



-( J F' 2 sin2F + 2F"sin 2 F)) =0. (II.1.26) 



It is independent of the dimension of the representation. Note that the differential 
equation is nonsingular only if F(0) — nir, n 6 Z. 

For the hedgehog form the baryon density reduces to the expression 



1 



3jV?r 2 

The corresponding baryon number is 

2 



j(.? + 1)(2j + 1) 



sin 2 F 



F'. 



B= / d 



,, Vn 3ti + !)( 2 i + ^(^(O) - 2 sin2F(0)). 



(II.1.27) 



(II.1.28) 



Combining the requirement that F(0) to be an integer multiple of n with the 
requirement that the lowest nonvanishing baryon number to be 1 gives the general 
expression for the normalization factor TV as 



N 



fj(j + l)(2j + l). 



(II.1.29) 



The equation of motion for chiral angle in the form (II. 1.26) depends on param- 
eter / w and e values. It is convenient to introduce a dimensionless variable f — ef n r 
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Classical 
chiral angle 



Figure II. 2. Classical chiral angle solution, as taken from Ref. 
(e = 5.45, = 129 MeV, f = eF^r). 



in which ( pl.1.26 ) takes the form 
F"(f)(l 



2Bin 2 F(f)^ +F/2( , ) sin2F(f) + 2 F/( , ) 



sin2F(f) sin2F(f)sin 2 F(f) 



0. (II.1.30) 



Numerical inv estig ation Q of (II. 1.30) leads to classical chiral angle solution F(f) 
shown in Fig. [L2, when boundary conditions F(0) = 7r, F(oo) = ensuring 



baryon number 1 are imposed. 

1.5. Higher order terms. Ther e exist s an infinite class of alternate stabiliz- 
ing terms for the Lagrangian density ( II.l.C ) , combinations of which can be used 
in place of Skyrme's quartic stabilizing term or be added to it |20|. An alternate 
term of quartic order (which for j = 1/2 yields the same result as the Skyrme term) 
is J2§ 



= -^5 (Tr{R M R„}Tr{R^} - (TrjR^}) 5 



(II.1.31) 



When this term is expressed in terms of the Euler angles (II. 1.3), the resulting 
Lagrangian density has the form (II. 1.17), with the exception that the stabilizing 
term that is proportional to e -2 has an additional factor + l)(2j + 1) |49|. 
Hence invariance of the physical predictions requires that the parameter 1 / e' 2 of the 
stabilizing term ( II. 1.31 ) be taken to be proportional to + 1) (2j + 1)) 

the parameter f n of the quadratic term to be proportional to + l)(2j + 1)) 
when a representation of dimension 2j + 1 is employed. Thus, Jz? 4 ' has different 



and 
l 



representation dependence from a similar term in (II. 1.6). 
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Consider then the sixth order stabilizing term |19, |2I 

J? 6 = e 6 Tr{[R li ,W}[R„ > TL x ][TL x ,TL IM }}. 
In terms of the Euler angles a. this Lagrangian density takes the form 



(11.1,32) 



_ j(j + l)(2j + l) 
-S6 — ^ e 



• 2 2 

sm a 



x d„a Zl dW 2 d v o? 3 d A a li d x a 15 d'V 6 . (11,1.33) 

This result reveals that the dependence on the dimension of the representation of 
this term is contained in the same overall factor j(j + l)(2j + 1) as in the Skyrme 
model Lagrangian (II. 1.17). Hence addition of the term J£§ maintains the simple 
overall dimension dependent factor of the original Skyrme model. 

As in the case of the qu artic ter m one can construct an alternative sixth order 
term, which is equivalent to ( II.1.32 ) in the case of the fundamental representation, 
but which differs in its dependence on j 



J ^v x v2v z Tr.fR R R T?' 7l R'? 2 'R ,?3 l 

^6 c /^?7l?72^3 - L1 \ Alj ^i- , - l '^2 /" 



(II.1.34) 



In terms of the Euler angles a this term also reduces to the expression ( |II.1.33[ ), 
with the exception of a n additi onal factor j(j + l)(2j + l)e' 6 /ee. Its dependence on 
j is thus different from ( II. 1.32 ) , although by adjusting the values of the parameters 
eg and ei differently in each representation equivalent dynamical predi ctions i n the 
classical fl model can be maintained. Obviously we cannot expres s ( II.1.31 ) and 
( [1. 1.34 ) as inner product of group generators, whereas for ( [1.1.33 ) this should be 
possible. 



2. The Lagrangian symmetries 

We start from construction of chirally invariant terms of the lowest order, which 
satisfy additional physical requirements. Maximal symmetry requirement together 
with B = 1 condition directly lead to Skyrme hedgehog solution in fundamental 
representation. This solution then is generalised to any SU(2) representation. 

2.1. Chiral symmetry breaking. It is commonly accepted that chiral sym- 
metry is the symmetry of QCD (theory of strong interactions) in the zero quark 
mass limit. There is only one second order chirally invariant term 

^(2)~ /^{R^R^d 3 *:, (II.2.1) 

whereas there are threej^ linearly independent invariants of the order four | 50 | 

JSf( 4 ) ~ a Tr{R M R M R^R"} + b Tr{R AI R l/ R^R l/ } + cTr{^R !/ ^R i/ }, (II.2.2) 

a, b, c being some constants. All three of them are Lorenz invariants so there are no 
reasons to prefer any one of them. But if we want to ensure positive energy density 
and (as a consequence) to avoid pathology in system dynamics we should take care 



9 This is not the case in quantum Skyrme model. 

10 Here we ana lyze only terms where all currents Rjj enter under single trace symbol. Gen- 
erally terms ( [1.2.2 ) and the term J Tr{R M R''} Tr{R„R"}d 3 x contribute to the same order in 
chiral perturbation theory fed]. 
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that time components of "right" currents entered the Lagrangian only in quadratic 
form. The requirement is satisfied by the only combination of the order fourpl 

if (4) = Tr{R Jli R' i R„R' y } - Tr{R jU R I ,R /1 R !/ } = ~ Tr{ [R„, R„] [R M , R"] }, 

(II.2.3) 



which is exactly the term suggested by T.H.R. Skyrme. The Lagrangian ( |LL6[ ) is 
invariant under global (point independent) chiral group SU(2)L eft ® SU(2)R ight of 
transformations of unitary field U 

U'(x) = VU(x)W- 1 ; V € SU(2) L ; W e SU(2) R . (II.2.4) 

The group SU(2)l®SU(2)r, however, is not a symmetry group of classical vacuum 
state (the highest symmetry field from Qo sector which takes on the constant value 
1) 

U(oc, t) = U„ ac (x, t) = 1, unless V = W. (II.2.5) 

As a consequence, maximal global invariant subgroup of configuration space of the 
model is 

diag(SC/(2) L ® SU(2) R ) » S[/(2) /sotopic , (II.2.6) 

where the standard notation diag(Gi ® G2) denotes the subgroup of Gi ® G2 with 
parameters of G\ and G2 being identified. 

2.2. Hedgehog ansatz as a lowest energy solution. In the baryon number 
zero sector the field which takes on the constant value 1 is the field of the highest 
symmetry. It is fully Poincare invariant and provides a classical description of the 
vacuum state. 

We expect that the ground state for B — 1 would be described by a configura- 
tion U with the maximal possible symmetry. When the winding number B is not 
zero, the field U cannot possess translational invariance. A translational invariant 
field is a constant and corresponds to B = 0. When Jj ^ 0, U cannot be rota- 
tionally invariant either. This is because a spherically symmetric unitary field U 
depends only on the radial distance r 

U+9 i U = x i Ut — , (II.2.7) 

where x is a unit vector in Cartesian coordinate system. Then 

#(U) ~ e ijk Tr{ (U+filfU) (U^U) (U+dfeU) } = 0. (II.2.8) 

To obtain one with B ^ let us blend isotopic rotations SU(2)i with space rota- 
tionsQ SO(3) to form a group diag(S77(2)i ® SO(3)). The unitary field U trans- 
forms under diag(5?7(2) I ® SO(3)) as follows Q: 

U ff (x) = T(p)U( ff - 1 x) = U(x); g G dmg(SU(2) I ® SO(3)). (II.2.9) 

Substituting expressions for the explicit rotation generators in the fundamental 
representation of SU (2)j yields differential equation 

-i(x x V).U(x) + [y,U(x)] = 0. (II.2.10) 



1X J?4 and JS?4 (II. 1.31) coincide up to the overall constant factor in the fundamental repre- 
sentation of SU(2) group. In SU(3) case these terms are different. 

12 SO(3) is homomorphic to SU(2). We keep notation SO(3) for spatial rotations (for a while) 
to make the separation more clear. 



20 



II. CLASSICAL SKYRME MODEL 



The solution^ of flH2.1C| ) is 

U(x) = U c (x) ee cosF(r) • 1 + ir ■ xsmF(r). 



(11,2.11) 



Pictorially it is ill ustra ted in Fig. [1,1 , where arrow length |sinF(r)| (for function 
F(r) itself see Fig. |ll.2| ) goes to zero as r — » oo an d r — * 0. The solution of ( H.2.10 ) 



is exactly the Skyrme's hedgehog ansatz ( |ll.l 2l|). G eneralization to the arbitrary 
representation is straightforward. Instead of ( H.2.10| ) we have 



i(x x V).U(x) + V2[Ji,U(x)] = 0, 



(IL2.12) 



where circular components (IF 1.8) are used both for the vector and isovector. The 
solution of (II.2.12) is a generalized hedgehog ansatz 



U c (x) = exp(-iv / 2 J a x a F{r)). (IF2.13) 

Sometimes the hedgehog ansatz ( [1.2.11 ), ( |IF2.13 ) is referred to as "spherically 
symmetric" solutions. These solutions are "spherically symmetric" only in the 
sense that a coordinate rotation is equivalent to an isospin rotation of the constant 
matrix A 

U c (x) = Aexp(-iV2 J a x a F(r))A\ A e SU(2). (IF2.14) 

To summarize, the highest symmetry solution for B ^ sector leads to certain 
mixing of indices associated with internal and geometric invariance (which are — a 
priori — completely unrelated). Similar examples are given by the monopole and 
instanton configurations which occur in SU(2) gauge theories. 

"One may wonder whether such a blend of internal and geometric symmetries 
may exist at a more fundamental level as a general feature of field theory and not 
simply in specific field configurations of particular models. This feature would be 
very attractive for the construction of a unified theory of all fundamental interac- 
tions including gravity. That this is not possible is expressed by the so called no-go 
theorems, in particular the theorem of Coleman and Mandula, which essentially 
says the following: the most general invariance group of a relativistic quantum field 
theory is a direct product of the Poincare group and an internal symmetry group, 
i.e. there is no mixture of these symmetry transformations. 

However, these no-go theorems do not claim that such a mixture cannot ex- 
ist if the set of all symmetry transformations represents a more general algebraic 
structure than a Fie group. Indeed, a famous result known as the theorem of Haag, 
Fopuszanski and Sohnius |87| states that the most general super Fie group of local 
field theory is the N- extended super Poincare group in which there is a non-trivial 
mixing of geometric transformations and internal SU(N) transformations. As a 
matter of fact, this result can also be viewed as a good argument in favour of the 
existence of supersymmetry as an invariance of nature since it states that super- 
symmetry is the natural (only possible) symmetry if one allows for super Fie groups 
as symmetry structures" Jg^ ]. 

The existence of nontrivial mixing for certain configurations in the Skyrme 
model as well as in supersymmetric models may serve as a strong argument for 
further investigations of the model which is much more simple to deal with than 
those of supersymmetric theories. 



The method of solution of equation (U.2.1C) is described, for example, in Ref. |5 
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2.3. Higher sectors solutions. It is proved |50[ (see, however, ||89|| ) that 
"spherically symmetric" hedgehog ansatz leads to the absolute energy minimum 
only when B = 1. For B = 2 value of baryon number it is hoped |50| that axial 
symmetric states 

F = F(r,tf); Q = 0(r,i?); $ = mv?; m G M, (II.2.15) 

realise energy minimum configurations. The statement was verified numerically (|90| , 
Bj]| . The value is Eb=2/Eb=i — 1-92 for the ratio of energies for axial symmetric 
solution of sector B = 2 and spherically symmetric hedgehog ansatz with B = 1. 
Energy /baryon densities for B = 2 configuration possess a toroidal symmetry (see 
Fig. |C.1| in Appendix ^|) . Stable ansatz for U minimizing energy and with baryon 
numbers B > 2 have been numerically found by various groups [17, 23, 92|. Energy 
densities for these static configurations have been plotted and a remarkable fact has 
been discovered that they are invariant under discrete subgroups Gr of the spatial 
rotation group SO (3). Some of them are shown [23] in Appendix^. The group Gr 
is the symmetry group of energy density. It is not necessarily the invariance group 
of the static U field. Published work [92 1 does not report on the symmetry group 
of the latter. 



Configurations with B > 1 are important in nuclear physics ||l7| since proton 
and nuclei could be related to quantized states of these soliton-like fields. Several 
recent studies support this point of view 



and suggest that the structures of 
heavier nuclei could resemble those of fullerene molecules, at least at the classical 
level. 

2.4. Reducible representations. Generalization of the model to arbitrary 
reducible representation is a bit straightforward. One needs only to sum over all 
irreducible repre sentat ions involving explicit dependence on representation. Thus, 
substitution for (IL1.4) is 



U(r,t) =^®D Jfc (a(r,t)). 



(II.2.16) 



The general scalar product ( II.1.9| ) then modifies to 

Tr(j-\J a Jb\j-) = I 1'"^^>U + l)(2j* + l)<J«,-6, 

k 

and the normalization factor ( |ll.l.29[ ) takes a form 

^ = |l>(ifc + i)(2j'fc + i)- 



(II.2.17) 



(II.2.18) 



Other formulas do not involve changes. 
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CHAPTER III 



Quantum Skyrme model 



This chapter contains main results. After brief remarks on quantization prob- 
lems in curved space we skip to collective coordinate approach and consider the 
Skyrme Lagrangian quantum mechanically ab initio. Assuming noncommutativity 
of dynamical variables we calculate expressions of Noether currents, magnetic mo- 
menta, axial coupling constant, etc. and numerically evaluate physical quantities 
using the classical chiral angle solution in various SU(2) representations. These 
numerical results then are used as starting input for self-consistent quantum chiral 
angle determination procedure. Numerical results of quantum chiral angle calcula- 
tions are presented in Appendices [X] and 



1. Quantization in curved space 

The purpose of the section is to remind readers the Dirac method of constrained 
quantization as well as problems of traditional quantization in curved space. The 
justification of the actual quantization method is considered without going into de- 
tails in the last subsection. The section, thus, provides the context for quantization 
procedure followed further but contains no new material. 

1.1. General remarks. Questions may be raised concerning the justification 
for quantizing the Skyrme Lagrangian at all, since it is not a fundamental field 
theory, but rather a classical model that results from taking the limit of such a 
theory, including only some degrees of freedom of the original theory. Nevertheless, 
there is a rich experience from the nonrelativistic many-body problems, for example, 



from nuclear physics [93], suggesting the validity of such an approach for the study 
of collective properties at low energies. 

The second remark concerns general quantization strategy. At the University 
level the construction of quantum theory passes three steps, namely Lagrangian 
(classical) — > Hamiltonian (classical) — » Hamiltonian (quantum) = Quan- 
tum theory. The quantization in collective coordinate approach [|4g|l , which we 
will follow in the work, slightly modifies this sequence. It starts from the quantum 
Lagrangian from the outset. By quantum Lagrangian we mean that dynamical 
coordinates q l and its time derivatives (velocities) q l do not commute. The explicit 
commutation relations at the moment are unknown. These relations are extracted 
from the standard commutation relations [q\p-') after we pass to quantum Hamil- 
tonian (and define canonical momenta p>). It can be shown that this modified 
formalism ||94) . 95 , 96] leads to consistent quantum description, which we will 



follow in subsequent sections. 

The results of modified and usual quantization sequences generally will dif- 
fer. Noncommuting quantum variables will generate additional terms while passing 
from quantum Lagrangian to quantum Hamiltonian. These terms are lost when we 
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impose canonical commutation relations after Hamiltonian is obtained. The prob- 
lem is similar to well known operator ordering problem. Further in the work we will 
refer to "canonical quantization method" when we start from quantum Lagrangian 
in the very beginning and to "semiclassical quantization method" when canonical 
commutation relations are imposed after Hamiltonian is calculated. It should be 



noted, however, that these two choices do not use up all possibilities. In Ref. |31|, 
for example, the following possibilities are discussed: 

-^Skyrme 

($,$) L hedg ,(F, F, A, A) £ co ii.(q,q) 

4 case I JJ- case II J) case III ■ (IH.1.1) 

-Hskyrme($,7r) => #hedg. {F, Uf , A, Jl A ) #coll.(q,p) 

In case I, the Hamiltonian gets by quantizing directly the Lagrangian in the Skyrme 
model from the beginning. After this quantization one can introduce the hedgehog 
ansatz and the collective coordinates. As another method (case II) we can do 
the quantization of the classical Lagrangian, which is obtained by introducing the 
hedgehog ansatz and after this we introduce the collective coordinates. Case III 
starts by getting the classical Lagrangian with the collective coordinates, using the 
hedgehog ansatz and getting the quantum Hamiltonian. The case II is free of the 
ordering problem among the operators. However, in both cases I and III, there 
are the ordering problems among the operators, and then quantum Hamiltonian 
cannot been determined uniquely. In summary, the problem is in which steps we 



do the quantization. Following ]48[ we utilize the following detailed quantization 
sequence: 



1. Introduce quantum collective coordinates A (q(t)) (III.2.1). They are quan 



turn in t he s ense that time differentiation of A requires the Weyl ordering 



(see Sec. |l,4|) . 



2. Make the Lagrangian quantum. Quantum dynamical variables are q(f) 



3. Following the method described in |48|, |9J, |95|, |96|, |97|, |98|, |99| pass to 

quantum Hamiltonian. 

4. Introduce the hedgehog ansatzQ, and solve integro-differential equation for 
quantum chiral angle. 

Another important point is symmetry properties of classical Lagrangian and 
quantum Hamiltonian derived from it. There exist quantization methods (for ex- 



ample, general covariant Hamiltonian method [(42|) preserving original classical La- 
grangian symmetries. The symmetric (Weyl) ordering of operators q and p (used 
in the work), however, cannot avoid a risk that the quantum Hamiltonian has chiral 



symmetry breaking term |31|. The feature thus can be used to provide (or explain) 
the origin of finite pion mass. 

The subsections below are intended to shed a bit light on the justification of 
the quantization method we will follow. 



Actually we introduce hedgehog ansatz before passing to explicit Hamiltonian. This is done 
only for reasons of simplicity of intermediate expressions and cannot affect the quantization itself. 
The only requirement of the quantization sequence in hand is that we do assume existence of 
solitonic solution ikql. 
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1.2. Towards quantum theory. There where attempts [ IOC ] to quantize 
classical equations of motion in the form 



n>^r = f, (III.1.2) 



dt 2 

without resort to a Lagrangian (or Hamiltonian) . The idea is to start with classical 



equations ( HI.1,2 ) and give commutation relations for the operators X 1 and X 3 

[X\X j ] = (III. 1.3a) 

[X\X j ] =iH6 i > j . (III. 1.3b) 



The result was that the existence of Lagrangian for (III. 1.2) essentially comes from 
( [II.1.3a| ). Thus Feynman's hope to quantize without a Lagrangian (or Hamil- 



tonian) w as do omed when he set the very reasonable condition that coordinates 



Ref. J102 1 . 



com mute ||10l|| . Let us briefly review traditional quantization methods following 



1.3. Traditional quantization methods. We will mention four of them here 
in the order of historical evolution. 

Heisenberg quantization is based on the following postulates (for the first pos- 
tulate see below): 

2. Introduce matrices Q = {Q mn } and P = {P mrl }, here m, n S {1,2,3, . . .}, 

that Satisfy [Q,P] mn = ^ (Qmp-Ppn PrapQpn) — ^^m.n* 

P 

3. Build a Hamiltonian matrix H = {H mn } as a function (e.g., polynomial) of 
the matrices, H mn = H(P, Q) m „, that is the same function as the classical 
Hamiltonian H{p 1 q). (In so doing there may be operator ordering ambigui- 
ties which this prescription cannot resolve; choose an ordering that leads to 
a Hermitian operator.) 

4. Introduce the equation of motion ihX mn — [X,H] m „ for the elements of a 
general matrix X = {X mn }. 

Along with these postulates comes the implicit task of solving the above equations of 
motion subject to suita ble op erator- valued boundary conditions. After the principal 



paper on quantization 103 . it subsequently became clear to W. Heisenberg that 



it is necessary to make this promotion from c-number to g-number varia bles only 



in Cartesian coordinates. Thus here is implicitly another postulate [ |104| : 

1. Express the classical kinematical variables p and q in Cartesian coordinates 
prior to promoting them to matrices {P mn } and {Q m n}, respectively. 



Schrodinger quantization assumes the following postulates 1 105 1 : 

1. Express the classical kinematical variables p and q in Cartesian coordinates. 

2. Promote the classical momentum p to the differential operator —ih(d/dq) 
and the classical coordinate q to the multiplication operator q, a choice that 
evidently satisfies the commutation relation [q, —ih{d/dq)] = iK. 

3. Define the Hamiltonian operator H as the classical Hamiltonian with the 
momentum variable p replaced by the operator —ih{d/dq) and the coor- 
dinate variable q replaced by the operator q. (In so doing there may be 
operator ordering ambiguities which this prescription cannot resolve; choose 
an ordering that leads to a Hermitian operator.) 
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4. For ipil) a complex, square integrable function of q, introduce the dynamical 
equation ifvtp = Hip. 

Implicit with these postulates is the instruction to solve the Schrodinger equation 
for a dense set of initial conditions and a large class of Hamiltonian operators. It is 
interesting to note that E. Schrodinger himself soon became aware of the fact that 
his procedure generally works only in Cartesian coordinates. 

Dirac quantization^ recipe deals with constrained dynamical systems. Constrains 
appear in the Hamilton formulation of all gauge theories we know of. Dirac- 
Bergmann constraint theory puts all these constraints into first or second classes. 
All second class constraints S a can be eliminated from the theory, whereas Dirac 
prescription for the implementation of first class constraints F a in quantum theory 
is that they be imposed as conditions on the physically allowed states |-): 

F a \-)=0. (III.1.4) 

Here F a is the quantum operator corresponding to the classical function F a . 
Feynman quantization focuses on the solution to the Schrodinger equation and 
postulates that the propagator, an integral kernel that maps the wave function 
(generally in the Schrodinger representation) at one time to the wave function at 



a later time, may be given by means of a path integral expression [106|. On the 
surface, it would seem that the (phase space) path integral, using only concepts from 
classical mechanics, could get around the need for Cartesian coordinates. That is 



not the case [107]. As postulates for a path integral quantization scheme we have: 

1. Express the classical kinematical variables p and q in Cartesian coordinates. 

2. Given that \q, t), where Q(t)\q, t) = q\q, t), denote sharp position eigenstates, 
write the transition matrix element in the form of a path integral as 

(q",T\q',0) = .# J exp((i/ft) / (pq - H(p, q)) dt) 9p 9q. (III.1.5) 



3. Recognize that the formal path integral (III.1.5) is effectively undefined and 



replace it by a regularized form of path integral, namely, 

. . N 

(q",T\q',0) = £m M N / exp(^(p I+1/a (<B+i - qi) 



c 



-H(p l+1/2 ,q l+1 + qi ))j UfL dp l+1/2 U^dqu (III.1.6) 

where q N +i = q", qo = q' , M N = {2i:h)~ {N+1 \ pi+ 1/2 = (pi +pi+i)/2, and 
e = T/(N + l). 

Implicit in the latter expression is a Weyl ordering choice to resolve any operator 
ordering ambiguities. Observe that the naive lattice formulation of the classical 
action leads to correct quantum mechanical results, generally speaking, only in 
Cartesian coordinates. Although the formal phase space path integral ( [II.1.5| ) 



appears superficially to be covariant under canonical coordinate transformations, 
it would be incorrect to conclude that was the case inasmuch as it would imply 
that the spectrum of diverse physical systems would be identical. In contrast, the 
naive lattice prescription applies only to Cartesian coordinates, the same family of 



2 The Dirac quantization although traditionally used is, actually, Hamilton formulation 
method, because operator ordering problems are not considered here. We give a very brief sum- 



mary of the theory in Sec. 
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[II. 1.6) of a path 



coordinates singled out in the first postulate of each of the previous quantization 
schemes. 

It is essential that traditional quantization methods start from global Cartesian 
coordinate systemP] and Hamiltonian function. But if we want Lorenz covariant 
theory, then it is not a good idea to start formulation from the Hamiltonian as 
there are no easy ways to ensure Lorenz covariance^. The easiest way to get Lorenz 
covariant theory is to start formulation from Lorenz invariant Lagrangian. This 
way leads us directly both to Hamilton formulation and operator ordering problems. 

1.4. Weyl ordering. A path integral formalism sometimes is referred to as 
a quantum method resolving the operator ordering problem [110 |. This is because 
the Weyl ordered expressions are used in the regularized form ( 
integral in order to resolve the ambiguities^. 

Let us illustrate the problem for an unconstrained system. Namely, let us recall 
the well known harmonic oscillator example. The problem is that two identical 
classical Hamiltonians 

H 1 ~x 2 +p 2 ; H 2 ~ (x- ip)(x + ip), (IIL1.9a) 

lead to different energies in quantum theory 

Ei ~ tuv(n + -); and E 2 ~ Hum, (HI. 1.9b) 

respectively, when the same commutation relations [p, q] is imposed. The reason 
is that the two quantum Hamiltonians ( [II.1.9a| ) differ exactly by the commutator 



[p,q], what is indicated by additional \frbj term in the system energy. The Weyl 
ordered Hamiltonian^ H\ results to the trueQ energy tuo(n + |). How can Weyl 
ordering can be applied in the general case? 



3 Needless to say that there are no ways to introduce global Cartesian coordinates onto arbi- 
trary configuration (phase) space. 

4 Lorenz covariant theory (Lorenz invariant Lagrangian) places time and space on the same 
footing, whereas time plays a special role in Hamiltonian and, therefore in quantum theory 



(quantum mechanics). Also, there are no ways to make time an operator [ LOfi [. Indeed, if time 
where an operator t it would be the component of a four-position operator X = (f, x) conjugate 
to the Hamiltonian H in the four momentum P = (H,p) 

[P»,X"]=ig» v ; ^" = {1,-1,-1,-1}. (III. 1.7) 

Then commutator [H, t ] = i implies 

cxp(- iei ) H cxp(iti) = H - e, (III.1.8) 

for any constant e. Thus the operator exp(iet) applied to any eigenstate \E) of H with energy 
eigenvalue E produces another eigenstate exp(iet)|_E) with shifted eigenvalue E — e. That indi- 
cates the presence of a continuous energy spectrum with range — oo < E < oo, contrary to the 
requiremen t th at E be bounded from below. Also, it contradicts the fact that generally, E is 



quantized [109]. 

5 Path integr al (IH.1.5| ) value still strongly depends on the finite-dimensional approximations. 



Points Pk,Qk m (HI. 1.6) are chosen in such a way (usually in the centre of the interval) that 



(HI. 1.6) limit coincides with operator eh tH matrix element only when i(pq + qp) corresponds to 



classical exp ressio n pg. Interesting, but the exceptions are known, when the middle point is not 
appropriate |107| |. As far as we know, there is no general recipe how to choose these points in the 
case of arbitrary curved space. 

6 The second Hamiltonian H2 is not Weyl ordered. Weyl ordering of it leads to the first one, 
because the Weyl ordering in two operator case is simple symmetrization. 

7 Although in this case the reference point can be shifted by — ^Hui to ensure the same result. 
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Because of nonvanishing commutator [q, q) ^ we need to state more explicitly 
what we understand under the symbol dtG(q). The most natural seems to be the 
definition 

U dGfeh 

« d e ti „ g 2 1 9 ' dq J 

which is a consequence of application of the Newton-Leibnitz rule to Taylor series 
expansion of arbitrary function G(q), 



(III.l.lO) 



G(q)=G(q Q ) + G'(q) q + Wiq) 

q=qo I 



q 2 + 



Indeed, 



(d t q )w = d t {qq) = qq + qq = -{q, 



1.. d(q 2 ) 



dq 



h 



(III.l.ll) 



(III.1.12a) 



3 1 d(n 3 ) 

(d t q 3 ) w = d t {qqq) = qq 2 + qqq + q 2 q = -(qq 2 + q 2 q) = -{q, -^i}, (IIL1.12b) 



(d t q n )w = d t (q . . q) = qq + ■ ■ ■ + qq = \ {(?, ^f-}. (III.1.12c) 

n n 

Here the notation (dtq n )w is related to the usual Weyl ordering notation (a n b)w 
in an obvious way 

^(d t q n ) w = (q n - 1 q) w . (III.1.13) 

The general Weyl ordered term (q n q)wi has a form 

1 " 

(q n q) Wl = —j E 9 n ~W- (HL1.14) 

It is straightforward to prove [ 1 1 1 1 that the above form of Weyl ordering is identical 
to the definition^ 

n | 

(q n q)w 2 = (-)" E w^ztk^' 1 ^ 1 - ( IIL1 - 15 ) 

The Weyl ordering has a number of interesting features [110] and is widely used. 
Further in the work we will follow definition (q n q)w 1 , which is identical to (III.l.lO) 
form. 

1.5. Hamilton formulation. Passage from Lagrangian to Hamiltonian some- 
times requires additional assumptions. Important class of theories, where standard 
Hamiltonization fails, is, for example, nonlinear models of elementary particles, 
including the Skyrme model. 

Let us concentrate on Lagrangian theories only, namely on local Lagrangian 
theories of the form 

L = y^'(0,9 AI 0)d™x, (III.1.16) 



8 Let us illustrate this for (q 2 q)w 1 and (q 2 q)w 2 terms. Indeed the sequence (q 2 q)w 2 
\{q 2 q + 2qqq + qq 2 ) = \ (? 2 <Z+(§ + Z)qqq+qq 2 ) = \{^q 2 q+^qqq+^qq 2 ) = \{q 2 q + qqq + qq 2 ) 
(q 2 q)w 1 shows the result. 
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allowing non-ambiguous equations of motionj]. All Lagrangians can be classified 
into two large groups, depending the system Hessian 

M, = (IH.1.17) 
is singular or not. In the case of nonsingular Hessian the usual method of Hamilton 



formulation is valid [112] 



1.6. Singular theories. When My is a singular matrix (detMy = 0), then 
we cannot express all velocities q as functions of momenta p and coordinates q. To 
show this we simply rewrite Lagrange second order equations 

SS B * ™* 0; S=[^t, (III.1.18) 



Sq l dq l dt dq l 

into two first order equations 

My v J = JQ , q l = v\ (III. 1.19) 

where the upper index v means that all q are changed to v. It is easy to see that the 
passage from Lagrangian to Hamiltonian is just variables exchange (q, v) — > (q, p). 
The exchange Jacobian is 

£>(q,v) _ 1 



D (l>P) det 



dv l dvi 



(III.1.20) 



and, therefore should differ from zero. It is a simple task to check wheather Skyrme 
Lagrangian is singular or n ot. T his depends on parametrization. In the Euler- 
Rodrigues parametrization ( II.l.l|) we can express one differential as a function of 



others, because of the constraint (|LL2|) on the fields. Thus we have singular theory 



in Euler-Rodrigues parametrization, whereas Hessian is clearly nonsingular in Euler 



angles parametrization ( II. 1.3 ) 



1.7. Dirac-Bergmann theory of constraints. |^| Let M be the space of 
"coordinates" appropriate to a Lagrangian L. We denote the points of M by 
q = (qi, 52, ■ ■ ■ )• Now given any manifold M, it is possible to associate two spaces 
TM and T*M to M. The space TM is called the tangent bundle over M. The 
coordinate of a point (q, q) of TM can be interpreted as a position and a velocity. 
The Lagrangian is a function on TM. The space T*M is called the cotangent 
bundle over M. The coordinate of a point (q, p) of T*M can be interpreted as a 
coordinate and a momentum. At each point q, momenta p belongs to the vector 
space dual to the vector space of velocities. Now given a Lagrangian L, there exists 
a map from TM to T*M defined by 



9q 

If this map is globally one to one and onto, the image of TM is T*M and we can 
express velocity as a function of position and momentum (see also previous section) . 
This is the case in elementary mechanics and leads to the familiar rules for the 
passage from Lagrangian to Hamiltonian mechanics. It may happen, however, that 

9 For example, Lagrangian ,£? = q leads to ambig uous Euler-Lagrange equation 1 = 0. 
10 The presentation of the subsection follows Ref. |ll3| . 
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the image oiTM under the map ([11,1.21) is not all of T*M . Suppose, for instance, 



that it is a submanifold of T*M defined by the equations 

Pa(q,p)=0; a = 1,2,.... (III.1.22) 

Then we are dealing with a theory with constraints. The constraints P a are said to 
be primary. 

The functions P a do not identically vanish on T*M: their zeros define a sub- 
manifold of T*M. A reflection of the fact that P a are not zero functions on T*M is 
that there exist functions G on T*M such that their Poisson brackets^ {G,P a } do 
not vanish on the surface P a = 0. These functions G generate canonical transfor- 
mations which take a point of the surface P a — out of this surface. It follows that 
it is incorrect to take Poisson brackets of arbitrary functions with both sides of the 



equations P a = and equate them. This fact is emphasized by rewriting ( [HI.1.22D , 
replacing the "strong" equality signs = of these equations by "weak" equality signs 
«: P a ~ 0. When P a (q, p) are weakly zero, we can in general set P a (q, p) equal to 

zero only after evaluating all Poisson brackets. 

In the presence of constraints, the Hamiltonian can be shown to be (TTij |Tl5| 

H = q a dL ^' q) - L(q, q) + V a P a {q, p) (III.1.23a) 

oq a 

= H + V a P a (q,p). (III.1.23b) 

In obtaining Ho from the first two terms of the first line, one can freely use the 
primary constrains. The functions V a are as yet undetermined Lagrange multipliers. 
Some of them may get determined later in the analysis while the remaining ones 
will continue to be unknown with even their time dependence arbitrary. 

Consistency of dynamic requires that the primary constrains are preserved in 
time. Thus we require that 

{P a ,H}^0. (III.1.24) 

These equations may determine some of the V a or they may hold identically when 
the constraints P a » are imposed. Yet another possibility is that they lead to the 
"secondary constraints" J^(q, p) » 0. The requirement {P^,i?} w may deter- 
mine more of the Lagrange multipliers, lead to tertiary constraints or be identically 



satisfied when ([II. 1.24) and P^~0 are imposed. We proceed in this fashion until 
no more new constraints are generated. 

Let us denote all the constraints one obtains in this way by Cj ~ 0. Dirac 
divides these constraints into the first and the second class constraints. First class 
constraints F a ~ are those for which {F a ,Cb} ~ 0, V6. In other words, the 
Poisson brackets of F a with C b vanish on the surface defined by Cb ~ 0. The 
remaining constraints S a are defined to be second class. It can be shown that 

{F a ,F b } =C c ab F c , (IIL1.25) 



where C^ h = — C§ a are functions on T*M. The proof can be found in JL13|, |114j1 . 
Let C be the submanifold of T*M defined by the constraints: 

C= {(q,p)\C b (q,Bp) = 0}. (III.1.26) 

Then since the canonical transformations generated by F a preserve the constraints, 
a point of C is mapped onto another point of G under the canonical transformations 



11 Recall that in quantum theory Poisson brackets are turned into commutators. 
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generated by F a . Since the canonical transformations generated by S a do not 
preserve the constraints, such is not the case for S a . Second class constraints can 
be eliminated by introducing the so-called Dirac brackets. They have the basic 
property that the Dirac bracket of S a with any function on T*M is weakly zero. 
Let F be the set of all functions which have zero Poisson brackets with S a - So long 
as we work with only such functions, we can use the constraints S a » as strong 
constraints S a — 0. Assuming that there are no first class constraints, the number n 
of functionally independent functions F is dim(T*M) — s, s being number of second 
class constraints. Thus s second class constraints eliminate s variables. Since matrix 
({S a , Sb}) is nonsingular and antisymmetric, s is even. Since dim(T*M) is even as 
well, n is even. 



Let us apply this theory to the Skyrme model |43|. In Euler-Rodrigues parametriza- 
tion (II. 1.1) we have primary constraint ( II. 1 .2) ) 



Si = $A - 1 = 0; a = 1,2,3,4. (IIL1.27a) 



The further requirement that the condition ( 111.1.27s] ) not vary in time can be 
satisfied by imposing a secondary constraint 

S 2 = $aM o6 $b - 1 = 0; a, 6 =1,2,3,4 , (III.1.27b) 

where M a b is inertia density matrix in Euler-Rodrigues parametrization explicit 



form of which is not important for further consideration (see |43] for details). It can 
be shown (using the canonical equations of motion) that this secondary constraint 
is independent of time [43]. These two constrains ( jIII.1.27 ) are of the second class, 



because Poisson bracket of the fields Si and S 2 is non-vanishing. This is also 
consistent with the fact that second class constraints come in pairs. 

If one carries out a canonical quantization by the usual commutators, trouble 
ensues from this noncommutativity [S , i,S f 2] pB 7^ 0. The quantum expression of 
the constraints is that every vector in Hilbert space must be annihilated by the 
constraint operators §1 and S 2 - It follows trivially that every vector must also 
be annihilated by their commutator and this conclusion is inconsistent since the 
commutator in question is itself nonvanishing in the canonical quantization. One 
resolution of this difficulty is to introduce modified classical brackets, the Dirac 
brackets, which share with the Poisson brackets all its basic algebraic properties, 
but are designed so that the Dirac bracket of any pair of second class constraints 
vanishes, in our case, [5 , i,S , 2 ] DB = 0. The Dirac brackets replace the Poisson 
brackets for determining the time evolution of relevant quantities. This replacement 
also eliminates the n eed t o introduce a Lagrange multiplier field. 

It can be shown 1 115| that for second class constraints there exists a canonical 



transformation leading to complete elimination of dependent dynamical variables. 
By canonical transformation here we mean non-singular transformation of dynam- 
ical variables 77 — > 77', 77 = (q, p); rj' = (q',p') if for arbitrary functions Gi{rf) 
and 62(77) Poisson brackets are invariant in the following sense 

{Gi(r ) ),G 2 (r,)} = {G' 1 (r,'),G 2 (r,')}, G'M) = Gifa), t = 1, 2. (IIL1.28) 

Note, however, that finding this eliminating transformation is nontrivial task, gen- 



eral solution of which is unknown |115[ . Dirac brackets in these new variables 



7]' coincide with Poisson brackets. In the SU(2) Skyrme model this elimination 
fortunately can easy be done by introducing Euler angles parametrization. 
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In SU(3) Skyrme model the problem cannot be completely solved by appro- 
priate parametrization, because of the first class constraints appearance in the 
theory. In this case Dirac prescription for dealing with first class constraints should 
be followed. Quantization of SU(3) Skyrme model is also complicated by Wess- 
Zumino term, which eliminates an extra discrete symmetry that is not a symmetry 
of QCD [29 , 30, 116]. Wess-Zumino term, however, vanishes in the SU(2) case 
even in quantum model |49|. 

To summarize, despite Dirac quantization^ (with its constraints classification 
into primary/secondary/tertiary/. . . , first/second class) is convenient and tradi- 
tionally often followed, it is, in fact not mandatory |112|, due to existence | 115| 



of canonical transformations (III. 1.28), which eliminate dependent dynamical vari- 
ables. If this elimination is technically formidable task, then Dirac procedure pro- 
vides us with consistent quantization method. Euler angles ( pi. 1.3 ) parametrization 
automatically eliminates dependent dynamical variables, and thus is consistent with 
Dirac prescription. 



2. Quantization of skyrmion in collective coordinate approach 

The section deals with quantization of the "zero frequency modes" or "collective 
coordinates" of the classical skyrmion. We proceed here (as in the classical case) 
with arbitrary irreducible representation. The case of reducible representation is 
investigated in the last subsection. 

2.1. Collective coordinate approach. An approximation of "zero modes" 
or "collective coordinates", which retains just a few modes out of a possible infinite 
number of modes, requires justification. It has been the subject of some criticism. 
We shall, however, proceed with our calculations using this approximation. 

Following G.S. Adkins et al. Q we shall employ collective rotational coordi- 
nates^] to separate the variables which depend on the time and spatial coordinates 

U(x,q(t)) = A(q(t))U c (x)At(q(t)), A(q(t)) G SU(2)j. (111.2.1) 

Three real independent parameters q(£) = (q 1 (t) , q 2 (t) , q 3 (t)) are dynamical quan- 
tum variables — skyrmion rotation (Euler) angles in the internal (isotopic) space 
SU{2) I = di&g(SU(2) L <giSU(2) R ) but not in geometric space SO(3). The skyrmion 
remains static in geometric space. 

Quantum fluctuations near the classical solutionis can be put into two different 
classes. Namely, fluctuation modes which are generated by action or Hamiltonian 
symmetries and modes orthogonal to the symmetric one. Symmetric fluctuation 
modes are of primary importance in quantum description because the infinitely 
small energy perturbation can lead to reasonable deviations from classical solu- 
tion. As a consequence, collective rotation matrices A(f), (t) in ( III.2.l| ) are not 
required to be small (i.e. close to the identity matrix). 



12 More precisely Hamiltonian formulation scheme. 

13 The method of collective coordinates originally was introduced by N.N Bogolyubov in |ll 
14 Recall that the most general collective rotational variables U(x, t) = V(t)U e (x) W _(t) 



corresponding to the full chiral invariance group SU(2)l ® SU(2)r of the Lagrangian (II. 1.6) 
are not appropriate, because we are interested only in fluctuations preserving U(x, t) = 1, when 
|x| — » oo. 
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2.2 . Commutation relations. We shall consider the Skyrme Lagrangian 
( [1,1.6 ) quantum mechanically ab initio. The generalized coordinates q(t) and 
velocities q(t) then satisfy the commutation relations^ 

[q r , q k ] = -ig rk (q). (III.2.2) 

Here the tensor <? rfc (q) is a function of generalized coordinates q only, the explicit 
form of which is determined after the quantization condition has been imposed^ 
The tensor g rk is symmetric with respect to interchange of the indices r and k as a 
consequence of the commutation relation [q r , q k ] = 0. Indeed, differentiation of the 
relation gives [q r , q k ] = [q k , q r ], from what it follows that g rk is symmetric. The 
commutation relation between a generalized velocity component q k and arbitrary 
function G(q) is given by 



[i\ G(q)] 



-V.9 fer (q)^G(q). 



(III.2.3) 



We shall employ the Weyl ordering for the noncommuting operators q, G(q) through- 
out 



(III.2.4) 



where symbol 9{G(q) further in the work is understood as (*9t G(q)) w . The operator 
ordering is fixed by the form of the classical Lagrangian ( |ILl,6j ) and no further 
ambiguity associated with it appears at the level of the Hamiltonian. In order to 
find the explicit form of g rk (q) one can substitute (III. 2.1) into (II. 1.6) and keep 
only terms quadratic in velocities! 17 ! 



L(q,q,F) =- /^(x,q(t),q(i),F(r))r 2 sinMrd^ 



(q)°-order term, 



N 

^\<TgrM)q k 



(III.2.5) 



where 



N=-j(j + l)(2j + l), 



(III.2.6) 



and we introduce here from the very beginning a normalization factor in the La- 
grangian ([1.1.6) in order to ensure baryon number 1. The 3x3 metric tensor 



<?r/c(q) is defined |49| as the scalar product of a set of functions (II. 1.14) Cy"^(q) 
orG; (m) (q) 



9rM) = -^(F)^(-) m G(™)(q)G(- m) (q) = -\a{F) ^(-)™G;M(q)G^(q) 



a(F)S r , 



1 

a(F)(8 rt i6k,3 + 5 r ,38k,i) cos g 2 , 



(III.2.7) 



15 We quantize only internal (isotopic) rotational degrees of freedom of the static soliton. 
Note that this does not imply that the quantization cannot affect the shape of the solution in 
geometr ic space. Conversely, because internal and geometric indices are mixed in th e solution 



(II. 2. 13) the shape of quan tum hedgehog ansatz is sig 

""Assumption (HI. 2. 2) actually is a consequence |94, 118| of canonical commutation relation 



( [II.2.1o| ) and canonical momentum definition ( jlll.2.9^ 



ca.ntly modified (see Fig. [II. 1) 



Lagrangian formulation and absence of the fir st order constrains suggest that terms linear 
in q should not appear. This is indeed the case | |48| . Terms which are independent of q do not 
contribute to momenta and thus to g rk also. 
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where a(F) (soliton inertia moment) is the following integral: 



a(F) 



d 3 x^(F(r)) = -3^-aOF) 



1 8tt 



dff 2 sin 2 .F(l + F 



sin 2 i^ 



*t. 3 i " V ' ' ' P (III2 - 8) 

The appropriate definition for the canonical momentum p r (which is conjugate to 
q r ) is 



Pr(q,q, F) 



dL(q,q,F) 1 , , ... 

- {9 ,flrfc(q)}- 



dg r 2 
The canonical commutation relation 

[p r (q.,q.,F),q k ] = -iS r k , 
then yields the following explicit form for functions g rk (q) 

ff rfe (q)=5 r - fe 1 (q). 
Note that for the time being we do not require [p r >J>k] = 0. 



(III.2.9) 

(III.2.10) 
(III.2.11) 



2.3. Angular momentum operators and remark on [PnPfc]- Because of 
the model spherical symmetry^ it is convenient to introduce operators J, J' instead 
of canonical momentum operators p 



J'n 



Ja 



{iV.C^q)} = (-1)'^ {? r ,C;(-«)(q)} , 
{jV,C[ a) (q)} = (-1)"^ {<f,C7<r<%)} . 



(III.2.12a) 



(m.2.12b) 



Straightforward but lenghty calculations then show that the following relations 
hold @: 

{pi,C( a) (<d},{p k ,C$ ) (q)}] = {^(^{C^q), [pi,Pk]}} 



2i 



1 1 



{p k ,C[ k a+b) (q)}. (III.2.13) 



.a 6 a + b. 

In obtaining ( ]III.2.13| ) the below equalities are useful: 

'{a, 6}, {c, d}] = {a, {c, [6, d]} } + {ft, {c, [a, d}}} + [a, {d, [6, c]}} 

+ {6,{d,[a,c]}}, (III.2.14) 

{a,6},c ={a, [6,c]| + {[a,c],&j, 



9 fc C; (o+fe) (q)-^C^ OJ (q) = 

[ft,G(q)] =-^G(q) 



V(o+ft)/ \ _ 



11 1 ■ 

a 6 a + &. 



C7^ ) (q)C; w (q), 



'(&)/ 



(III.2.15) 

(III.2.16) 
(III.2.17) 



18 Important here is the spherical symmetry of dynamical degrees of freedom q : SU (2)j 
SO(3)/. 
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where G(q) is an arbitrary function of q. Thus from ( [11,2.13 ) we have 



[%, J' b ] = - j({C('i)(q), {<$)(<!), [pi.Pft] }}^ 





1 


1 1 ■ 


)+ 


.0 


b a + b. 



J' b . (III.2.18) 



Let us examine commutator [p r >i>fc]- From ([II. 2. 12a) we have 

Pr = £{Cr (0) (q),£}- (HI.2.19) 

Introducing vielbeins for (q) and dual vielbeins /i^s (q) 

^(-l)^W(q)4-°)(q) = ffrfc (q), ^(-l)^[ a) (q)/if_ a) (q) = ^(q), 

(IU.2.20a) 

4 a) (q)=^y^^(q), 



^(a)(q) = -* 



a(F) 



Cfo(q), (III.2.20b) 



we can rewrite ( [II.2.1S ) in a geometrically more suitable form 



Pr 



Utilizing equations 



[j>< o) (q)],^(q) 
[J^4 a) (q)],[j^^)(q)] 



= o, 

= 0, 



(III.2.21) 

(III.2.22a) 
(III.2.22b) 



which are obtained without reco urse to |p r ,pfc] we obtain with the help of (HI. 2. 21), 
( EII.2.14 ), CTII.2.16D , ( |III.2.15D , ( |II.1.15b| ) and ( ]III.2.17| ) 
1 



[Pr,Pk] = 



2a(F) 



{^ a >(q),ja.{4 6) (q).^} 



^(^{e)(q),{4 b) (q),[j:,^]}}+2{j^[j:,/ l ( b )(q)]4 a) (q) 

-[J^e)(q)]4 b) (q)} 

^{4 a Hq),{4 b) (qM«]}} 



2a(F) 
1 

+ 4 



I J' 



11 1 

m n m + n 



C7;( m )(q)C7f l) (q)}. 



(III.2.23) 



I t is easy to see that metric flHI.2.7j ) is invariant under local rotations of vielbeins 
( |EII.2.2cj ) 



/ l W(q)^/ l ;W(q)= J D^(a(q))/ l W(q). 



(III.2.24) 



This is related to the fact that space defined by metric (HI. 2. 7) is of constant 
curvature 

3 



R = R r 



kr 



2a(F) ' 



h 



sh pi 



T r , h T 



lh ps J 5 



(III.2.25a) 
(HI.2.25b) 
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r r , being Christoffel symbols 



rr u = ^9 rp (dig P s + d s g pi - d p gis 



(111.2.25c) 



Equatio ns ( III. 2.2 ) and ( IIL2.ll] ) then imply that one can locally rotate vielbeins 
( [II. 2.24 ) without affecting commutation relation ( HI. 2.10 ). 

Defining spin connection A a b . r (q) in usual way [4S, 97, 119 ] 

V r C' k (a } (q) =C^ (a) (q) - r s rfe C^(q), 



v4 o) (q)=-iV(q)cf(q), 



(III.2.26a) 
(III.2.26b) 

and employing (111.2.25c) and ( [III. 2.16 ) we get explicit form for connection A a b . r (q) 

1 



i i r 

b a — b a. 



c; (b) (q), 



(IIL2.27) 



which transformes under local rotation of vielbeins ( |III.2.24| ) in well known way 

A\M) = ^, c (a(q))A c d . r (q)£>^( a(q) ) _ ^£>l jC ( a (q))l>*t t (a(q)) . (IIL2.28) 
Rewriting equation ( |IIL2.23| ) in terms of spin connection ^4 a b . r (q) 

or in a form reflecting the expression symmetry in respect of interchange of indices 

r, k 



2^){^ )(q) '^^ )(q) '[ J - J ^} 



[j' m , A m „. r (q)4">( q ) - A"\ l;fc (q)/4")(q)}, 



(III.2.29) 



we are ready to show that one can perform a suitable local rotation ( [HL2.24P of 
vielbeins around a point P(q) to ensure [p r ,pfe] = at that point. From Jacobi 
identity 



we have 



[£. [v,C]] + [v, [££]] + [C, [Z,ri\]=o, 



[[Pr,Pk],G(q)] = 0, 



(III.2.30) 



(IIL2.31) 



for an arbitrary function G(q) satisfyi ng d r dkG(q) = dkd r G(q). Assume that 
such a rotation has been found. Then ([II.2.1g|) reduces to 



[J'a > j'b\ — 



1 1 



1 



a b a + b. 



J a+bi 



(III.2.32) 



and ( III.2.23 ) to 



\Pr,Pk 



[j' m , A m „ ;r (q)4")(q) - A^q^q)}, (III.2.33) 
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wher e we hav e used symmetry reflecting form ( III,2.29| ), Substitution of ( QI.2,33 ) 
into (|HI.2.31|) leads to 



(A ro n . r (q)^(q3 - A m n;k (q)hi n \q))hl m) d s G(q) = 0. 



(III.2.34) 



This implies that A m n . r (q) ~ in the very vicinity of the point -P(q). Thus, 
equation ( [11,2.28 ) with A' a b . r (q) = gives us explicit partial differential equations 
for local rotation angles a(q). 

When [priPfcl = 0, opera tors J', J become angular momentum operators with 
usual commutation relations ( ITII.2.32D and 



[j a , Jb] 



1 1 

a b 



1 

a + b 



(III.2.35) 



respectively. The operator J' is then a "right rotation" generating matrix D (q) 



[4i ^m,m'(q) 



= - £ 



Jn 



and J is a "left rotation" generating matrix D £ (q) 



■'"I -^m.m'(q) 



I , m 



Ja 



m-a)D l m _„ m /(q). 



(III.2.36) 



(III.2.37) 



2.4. The Hamiltonian and state vectors. In order to calculate all terms 
in quantum Lagrangian expression ([II. 2.5), we substitute collective rotational co- 
ordinates (III. 2.1) into properly normalized Lagrangian density ([1. 1.6). Utilizing 
commutation rule ( III. 2.3] ) and relations ([1. 1.12) for Wigner D J functions we can 
pull out all velocities to one or another side in symmetric fashion. Formulas^ 



Tr{(jm\j' a Ji\jm)} = (-l) a -j(j + l)(2j + l)J a ,_ 6) 



Tr{(jm\j' a J{,J' c \jm)} = -(-1) 



a j(j + l)(2j + l) 
3-4 



1 1 

c 6 



(III.2.38a) 
(III.2.38b) 



Tr{(im|J^44|im>} = ^(-ir+ 6 ^ 2 (j + l) 2 (2i + l) 2 {^ { 

k 



l l 



k 

a + b. 



11 k - 

c d c + d. 



(III.2.38c) 



then allow us to take trace explicitly. The resulting expression still contains a lot 
(up to 6 for Lagrangian density and up to 8 for Noether current densities) of sums 
over repeated group indices. We have used computer algebra system to make the 
explicit summation^. The result of all computations is the following explicit form 
for the quantum Lagrangian j2|: 

L(q, q, F) = —M(F) - AMj(F) + J' 2 , (III.2.39a) 



-M(F) — AMj(F) 



a(F) 
1 

a(F)' 



(III.2.39b) 



19 Formulas ( IH.2.38| ) are derived by E. Norvaisas (private communication). 

20 Recently we have made some progress in completing the summation manually, but inter- 
mediate results are still very large. Computer algebra system also was used to check many of 
symbolic manipulations mentioned above. 
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and its density 

^q ) q) = ^(j' 2 -(J'-x)(J'-x 
Here M(F) is classical soliton mass 

e 



A^j(F) - Jt(F). (III.2.40) 



M (F) = 



d 3 xJ?(F(r)) = J —M(F) 



-Ztt 



dff 2 I F 



sin 2 F 



2F' 



sin 2 F 



(III.2.41) 



and 

AM^F) 



d 3 xA.4(f(r)) = e 3 /iAMj(F) 



2ne 3 U 



5o : 



r,2 



J dff 2 sin 2 F^5 + 2{2j - l)(2j + 3) sin 2 ^ + (8j{j 



+ (2j(j + 1) + 1) 



sin 2 .F 



- 2(2j - l)(2j + 3)F' 2 sin 2 F 



1) - l)F 



(III.2.42) 



i s quant um mass correction. The angular momentum operator on the r.h.s. of 
( PI.2.4CD can be separated into scalar and tensor terms in the usual way 



(J'.x)(J'-x) 



47T, 



-Yr 



2,m+m' & <f) 



1 1 2 

m m! m + m' 



V V , 

(III.2.43) 



where l/, m (i?, ip) |82| is a spherical functions. 

The volume integral of the Lagrangian density ( |[II.2,4C|) gives the Lagrangian 
( [II.2.3S|). In th e fundamental representation (j — 1/2), the second rank tensor 
part of JlII.2.43j ) vanishes. This implies that the quadrupole moment of the A 33 
resonance canno t be described within the fundamental representation. 

It is known |ll8(| that in the quantum mechanics, the Hamilton formalism is 
inconsistent (see also Sec. |]) with Lagrange one for velocity dependent potentials 
if the Hamiltonian is defined by the ordinary method 

hp r ,f}-L, (III.2.44) 



K 

and that this K does not satisfy the canonical equations of motion 
dH(q,p) . t dH(q,p) 



q ; 



(III.2.45) 



with H replaced by K. In a number of works Q |95|, |9^, |08|, |S9|, |l!8| , |l20| , fm], 

122 1 it has been shown that co nsistent Lagrange and Hamilton formalism exists^] 
for constant curvature ( [II.2.25| ) spaces. The correct Hamiltonian is given ||48]| by 

1 



H = -{p r ,q r }-L(c l ,q)-Z(c l ) 



(III.2.46) 



satisfying the canonical equations ( [III.2.45 ) of motion. The extra term Z(q) arises 
from noncommutativity of operators and can be expressed in terms of g rk (q) and 
g r fc(q), so that it does not involve q or p. There exist a few explicit forms of Z(q). 

21 Moreover, consistent variation in connection with Noether theorem exists for the case of 
constant curvature l|98l. 
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Calculations become shorter if one uses the following expression of Z(q), given by 
R. Sugano 



^9 rk (dr9 sp )(d k g sp ) + l -g rk (d r g^)(d s9kp ). 



By direct calculations we see the following identity to hold: 



Thus the true Hamiltonian is written as 

H 3 (F) —M(F) + AMj(F) 

—M(F) + AMj(F) 



I'2 



a(F) 



a(F) 



(III.2.47) 

(III.2.48) 

(ni.2.49a) 
(III.2.49b) 



Note, that both the second and the third terms on the right-hand side of ( [HI.2.49D 
are of the order of h 2 , the magnitude b eing cha racteristic of operator ordering 
contribution [35, 42]. The Hamiltonian ( III. 2.49 ) yields |^^| canonical equations 

of motion ( [II.2.45 ) under the condition J = J' = 0. This condition is consistent 
with H (q, p) to be a Hamiltonian, because we have 



[J',#(q,p)] =0; [J,ff(q,p)] = 0, 



(III.2.50) 



from commutation relations (III. 2. 32) and (111.2.35). The most important feature 
of this result is that the quantum correction AMj (F) is negative definite and that it 
depends explicitly on the dimension of the representation of the SU (2) group. This 
term is lost in the usual semiclassical treatment [Q of the Skyrme model even in the 
fundamental representation [48 1 of SU(2), because that ignores the commutation 
relations. For the Hamiltonian (III.2.4S) the normalized state vectors with fixed 
spin and isospin £ are 



Wisospin, m sp \ n 

These have the eigenvalues 



V27TT 

47T 



H(j ) £ ) F) = M(F)+AM j (F) 



1(1 + 1) 
2a(F) 



(III.2.51) 



(IIL2.52) 



This expression is the quantum version of the mass formula of the Sk yrme model. 
The Hamiltonian density corresponding to the Hamiltonian ( III. 2.49 ) has the fol- 
lowing matrix elements for baryon states with spin and isospin I > 1/2: 



I 

mi, m s 



^(r,q) 



t 

mi,m s 



=J?(F{r)) 



(F(r)) 



2a 2 (F) 

^(3m 2 s -e(l + l))Y 2 , {#,p) 



l+l) 
(III.2.53) 



For nucleons I = 1/2 the dependence on angles is absent and the quantum skyrmion 
is, therefore, spherically symmetric as required. The states with larger spin than 
1/2 are, thus, described as deformed. 
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2.5. Reducible representations. Here we provide modifications which are 
necessary when passing to the general reducible representation 

D 3 '->5^ffiD A . (III.2.54) 



Thus AM j (|[II.2.42D in quantum Lagr angian (HI.2.39| ) (and its density (|HI.2.40D ) , 
Hamiltonian ( p!II.2.49|) (and its density ( [11.2,53 )) and quantum mass formula ( III.2.5S ) 



is modified to 

AM^(F) = J d 3 xAJ' Ej (F(r)) = e 3 /,AM s .(ii') 

= - 2 -^fl d ^ 2 sin 2 ^(l5 + 4d 2 sin 2 F(l - F' 2 ) + 2d 3 ^- + 2d 1 F' 2 

(III.2.55) 

The coefficients di,d 2 , d 3 in these expressions are given as 

di = ^^3k{jk + l)(2j fe + l){8 Jk {j k + i) - i), (III.2.56a) 

fc 

d 2 = ^^jfcO'fc + l)(2jfc + l)(2 Jk - l)(2 Jk + 3), (III.2.56b) 

fc 

d 3 = ^^ife(ife + l)(2i* + l)(2j fc (7* + i) + 1), (111.2.56c) 

fc 

with the same generalized normalization factor as in classical model ( 11.2.1^ ) 

N = Ij^JkUk + l)(2i fc + 1). (IIL2.57) 

d fc 

Further in the work we proceed with the general reducible representation. 

3. The Noether currents 

In the section we introduce vector and axial-vector transformations of unitary 
field U and calculate vector and axial-vector currents assuming non-commutativity 
of dynamical variables in the Lagrangian from the outset. Obtained expressions 
then are used to define nucleon and A33-resonance magnetic momenta and axial 
coupling constant $a- We show that in this approach axial symmetry becomes bro- 
ken, whereas vector symmetry is still conserved. In the end we provide definitions 
of various experimentally measurable radii following Ref. |^| . 

3.1. Vector and axial transformations. Long before QCD was believed 
to be the theory of strong interactions, the phenomenological indications for the 
existence of chiral symmetry came from the study of the nuclear beta decay: n — » 
p + e + D. There one finds that the weak coupling constants for the vector and axial- 
vector hadronic currents, V and A, did not (in the case of V") or only slightly(25% in 
the case of A) differ from those for the leptonic counterparts. Consequently, strong 
interaction "radiative" corrections to the weak vector and axial vector "charge" are 
absent. The same is true for the more familiar case of the electric charge and there 
we know that it is its conservation, which protects it from radiative corrections. 
Analogously, we expect the weak vector and axial vector charge, or more generally, 



22 The presentation of the subsection follows Ref. |85 
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currents, to be conserved due to some symmetry of the strong interaction. In case 
of the vector current, the underlying symmetry is the well known isospin symmetry 
of the strong interactions and thus the hadronic vector current is identified with the 
isospin current. Vector-transformation, therefore, is the isospin rotation. In terms 
of pions this can be written as 

7r^7r + 0X7r, (III.3.1) 

which states that the isospin direction of the pion is rotated by angle 0. The 
isospin, thus, is a constant of motion associated with vector transformation. Also 
note that vector transformation leaves the vacuum U = 1 invariant, due to rotation 
in the isospace only. 

The identification of the axial current, on the other hand, is not so straightfor- 
ward. This is due to another, very important and interesting feature of the strong 
interaction, namely that the symmetry associated with the conserved axial-vector 
current is "spontaneously broken". By that, one means that while the Hamiltonian 
possesses the symmetry, its ground state does not. An important consequence of 
spontaneous breakdown of a symmetry is the existence of a massless mode, the 
so-called Goldstone boson. In our case, the Goldstone boson is the pion. If chiral 
symmetry were a perfect symmetry of QCD, the pion should be massless. Since 
chiral symmetry is only approximate, we expect the pion to have a finite but small 
(compared to all other hadrons) mass. This is indeed the case. Interpretation of 
axial-transformation is also not so straightforward. It can be shown [|8j| that this 
transformation mixes pion and a meson states 

7T — > 7T + 0(7, 

a — > a — • 7r. 

The pion thus is "rotated" into sigma meson under the axial transformation and 
vice versa. 

Let us consider the matrix element of the axial current between the vacuum 
and the pion: (O^^tt). Because of parity, the matrix element describes the weak 
decay of the pion and must be proportional to the pion momentum (this is the only 
vector around) 

<0|^(x)|7r b (p)) = -iUp^6 a b e-^ x , (III.3.2) 

where p is pion momentum, indices a, b refer to isospin, \x indicates the Lorenz 
vector character of axial current, and f n is the pion decay constant determined 



from the experiment^. Taking divergence of (III. 3. 2) we obtain the relation which 
is often in the literature referred to as the PCAC relation (partial conservation of 
axial current) 

<0|9^(x)|7r fc (p)) = -if^mlS^e-^ = -^.m^^"''. (III.3.3) 

From this equation we see that to the extent that the pion mass is small compared to 
hadronic scales, the axial current is approximately conserved. Or in other words, 
the smallness of the pion mass is directly related to the partial conservation of 
the axial current, i.e. to the fact that the axial transformation is an approximate 
symmetry of QCD. 



23 Instead of e lp x generally one should use / p (x) — the asymptotic pion wave function as 
a solution of the Klein-Gordon equation ]35||. 
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3.2. Vector and axial currents in the Skyrme model. The Lagrangian 
density of the Skyrme model is invariant under left and right transformations of 
the unitary field U(q) 



U -> (left)U = (1 - i2y/2u a J' a )V, 
U -> (right)U = U(l + i2V2uj a J , a ). 



(III.3.4a) 
(IU.3.4b) 



The vector and axial Noether currents are nevertheless simpler and directly related 
to physical observables. They are associated with the transformations [0] 



U(x) 



Vector transf. 
Axial transf. 



(l - i2V2uj a J'^j U(x) (l±i2V2u a j' a ^ , (III.3.5) 



respectively. T he fact or —2^/2 before the generators is introduced so that the 
transformation (111.3.5) for j = 1/2 matches the infinitesimal transformation in 
Ref. The corresponding Noether currents can be expressed in terms of the 
collective coordinates (III. 2.1). After this substitution Noether currents become 
operators in terms of the generalized collective coordinates q and the generalized 
angular momentum operators J' ( 111.2.12a ). Long manipulations similar to those, 
described in obtaining quantum Hamiltonian, lead to the explicit expression for the 
vector current density 



-a d^fy _ 

2v'2sin 2 F 



(left) 



d ( W>) 
1 



(right) 



sin 2 F 2d 2 + 5 



sin 2 F 



sin 2 F 
V2e 2 a 2 



r 2 4 • 5 • a 2 
(-!)■( [3 1 vx]_X s (q)[[j'xx]xx 

[j'xx] xxl X s (q)[j'xx], 



111 

u s b 



(III.3.6) 



Here V b is a circular component of the gradient operator. The indices a and b 
denote isospin and spin components respectively. The time (charge) component of 
the vector current density becomes S 



* d(dauj a ) 



(left) 



2V2(-iy 



8 (y UJa) 

1 



+ (right) 



d(V UJ a 

sin 2 F 



(q)z s (J'-x 



(111.3,7) 
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The explicit expression for the axial current density takes the form 



b ~d{V b u a ) 



suv2F 



(left) 



as? 



=[fi 



d(V b UJa 

1 sin2.F 



(right) 



d{V b LUa) 



(f 12 



&\tl 2 F sin 2 F 



4a 2 



Dl b (q)+(f 2 (2F> 



sin2F\ 



1 /F /2 sin2F sm 2 Fsm2F AF'sm 2 F sm 2 Fsin2F 



e 2 ( 



e" \ r 
IF 1 sin 2 F(-l) s 



Aa 2 r 



DlMx.J' 2 + 3 !2 DlM)x. s - 2D* (q)z_ s (J' ■ x)(J' • x) 



sin ^sir^F 



(-lW [j'xx] xx D^(q) [j'xx] xx 



[J' x x] x x Dl s (q) [J' x x] x x 



(III.3i 



The operators ( pi.3.6[ ) , ( |III,3.7| ) and ( |III.3.q ) are well defined for all representations 
j of the classical soliton and for fixed spin and isospin £ of the quantum skyrmion. 
The new terms which are absent in the semiclassical case are those that have the 

factor a (F) in the denominator. 

The matrix element of the divergence of the vector current density ( HI.3.6| ) 
vanishespl 



£ 

uii,m s 



rrii, rrii 



= 0. 



(III.3.9) 



The re sult just confirms validity of the variation procedure on constant curvature 
space 1 1 20 1 . 



3.3. Baryon current density and magnetic momenta operators. The 

conserved topological current density in the Skyrme model is the baryon current 
density. For the hedgehog solution its components take the form Q 

&Jx,F(r)) = —= F'sin 2 -F[J' x x] . (III.3.10) 

It is sketched in Fig. frll.l for classical and quantum chiral angles. 

Rot ating sol iton generates isoscalar magnetic moment associated with current 
density ( III.3.1C ). The matrix elements of the third component of the isoscalar 
magnetic momentum operator have the form H 



i 

rrii, rn s 



£ 

rrii, m s 



£ 

mi,m s 



d 3 xr [x x £S] 



£ 

rrii,m s 



,1/2 



e {£(£ + !)) '- . 2 ,r£ 1 £ 



(III.3.11) 



■in 



m s m 5 



24 The v ector-transformation (HI. 3. 5) (finite) U — » BUB -1 and collective coordinate rota- 
tions (HI. 2.1) possesses the same symmetry, whereas axial-vector transformation (finite) U — > 
BUB does not. 




Figure III. 1 . Baryon charge density distribution. Dashed line 
denotes classical chiral angle distribution, with parameters e = 
5.45; /„■ = 64.5 MeV taken from Ref. J^|. Quantum chiral angle 
distributions for various j values are plotted with solid lines. 



where the isoscalar electric mean square radius is given as 



(rh=o) 



2 fi 



,j=o); (?e,i=o) — 



r 2 



F'(r)sin 2 F(r)dr. (IIL3.12) 



The matrix elements of the third component of the isovector part of magnetic 
momentum operator that is obtained from the isovector current ( HI. 3.7 ) have the 
form 



I 

mi,m s 



[W=i]j 



t 

mi,m s 



~ / d 3 xr[x 



a— 3-, 



a e an 
e 3 U /tt 3a 2 



df f 2 sin 4 F ( 1 



I 

morris 

1(1+1) _ _ch_ 
3 2-5 



+ 



(-1) 2<? ( M(l+l){2l-l)(2t+\)(2l+S)\ a f 1 2 1 



2-:! 



)'{ e }) 



- £ 1 


£ - 




■ £ 


1 


£ ■ 


m s 


m s _ 




mi 





mi. 



(III.3.13) 



where the symbol in the curly brackets is a 6j coefficient | 82 | 
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From ( III. 3.11 ) and ( III.3,13[ ) proton and neutron magnetic moments measured 
in nuclear magnetons can be extracted using relations 



Pp —2 ((l/2a%l^ 7 = lHl/2a/2) + (l/ 1 2(l 2 /2|[Al=l]3| 1/ 1 2 ^ 1 2 /2 )^ , 
Mn =^ ((-l/2 2 l/2|[/'J = 0]3|-l/2 2 l/2)-(-l/2 2 l/2|[Al=l]3|_ 1 /2 2 l/2) 

Similar formulas exist for A33-resonance magnetic moments 
1 



A 



33 



A 



2 
1 

: 2 
1 

2 
1 

33 =2 



33 



A 
^33 



,3/2(3/2! [Ar=o] 3 1 3/2(3/2) + (3/2(3/2 1 [Ar=i]s 13/2(3/2) 

2 ((1/2(3/2! [Al = o]3 1 1/2(3/2) + ( 1/2(3/2! [A/ = l]3 I 1/2(3/2)) 1 



t - 1/2, 3/2 1 [£/ = <>] 3 I -1/2.3/2) + (-1/2,3/2! [A/=l]s 1-1/2,3/2/ 



3/2. 3/2 1 [Ar=o]3 I -3 /0.3/2) + ( 



-3/2 2 3/2|[A'=l] 



3/2 \ 
J I -3/2,3/2/ 



(m.3.14a) 
(III.3.14b) 

(in.3.15a) 
(IIL3.15b) 
(111.3.15c) 
(IIL3.15d) 



Relations ( III.3.14 ) also allow us to obtain proton and neutron charge distribu- 
tions, when isoscalar and isovector current densities^] are used instead of integrated 
expressions (III. 3. 11) and (III. 3. 13). These distributions are measurable quantities 
for which semi-empirical formulas exist [ 123] 

1\ J3 ^2 



Pp = 



■ exp(-M D r) ; M D = 0.84 GeV, 



(III.3.16a) 



~lJL n M%r ( 8m 2 N (e-Kp(-M D r) - exp(-2m N r / \fb£)) 



2(5.6M|, - Am%) \ 



5.6M|, 



Am 2 N 



tMd exp(— M^r) 



(III.3.16b) 



Distributions for classical (short-dashed) and quantum chiral angles in v arious rep- 
resen tations (s olid lines) are plotted vers us se mi-empirical distributi on (III.3.16b ) 
and ( 111.3.16a ) (long-dashed line) in Fig. III. 2| for proton and in Fig. QI.3 for neu- 
tron. The neutron distribution peak value (<~ 0.5) in classical Skyrme model is 
known to be too large when compared with the empirical value (0.2 — 0.3). Our 
results for quantum chiral angle show significant improvement in neutron charge 
density distribution, which only weakly depends on the representation used. Note 
also that measurement indicates much faster distributi on fa ll in the asymptotic re- 
gion even when compared with semi-empirical formula 1 123 . Speed up of quantum 
chiral soliton distribution fall t hus is also a step in right direction. For experimental 
distributions we refer to ||l23|| and references therein. The axial coupling qa mea- 
sures the spin-isospin correlation in the nucleon and is defined as the expectation 



value of the axial current A a in a nucleon state at zero moment um tran sfer |53 |. 
Comparison of the Fourier transform of the axial current density ( III. 3.8 ) with the 
usual phenomenological expression for the axial-vector current leads to the ax ial- 
vector coupling constant expression^]. From the axial current density ( HI. 3.8 ) we 
obtain the axial coupling constant qa of the nucleon as H] 



9A 



1/2 
1/2,1/2 



1/2 
-1/2,1/2 



^ 2 e 2 
3a 2 (F) 



(III.3.17) 



25 These densities are integrated over angular variables tp, i9. The factor r 2 (which comes from 
the Jacobian) is also included to ensure usual dimensions. 

26 In the limit of vanishing pion mass, the resulting expression for gA should be multiplied 
by a factor pi |. There are subtle points in calculating this quantity for solutions with ~ 1/r 2 
asym ptotic g|. This factor should not be used in calculating gA in finite pion mass model |8|, p5| , 
123], and thus in quantum case either. 
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Figure III. 2. Proton charge density distribution. Short-dashed 
line denotes classical chiral angle distribution, with parameters 
e = 5.45; = 64.5 MeV from Ref. Quantum chiral an- 
gle distributions for various j values are plotted with solid lines. 
Long-dashed line denotes semi-empirical charge density distribu- 
tion ||123||. 



where 

gx (F) = y J df [r 2 F' + r sin 2F + rF' 2 sin 2F + 2F' sin 2 F + ^7-^ sin 2f) . 

(III.3.18) 

All versions of the topological soliton model lead to underprediction of gA com- 
pared to the empirical value of 1.26. In the simple Skyrme model with a pion mass 
term | 132| | the predicted value for gA was only 0.65. By readjusting the param- 



eters so as to fit empirical values for / w and ttNN coupling constant the value 
is somewhat increased (to 0.82). In the vector-meson-stabilised model |124| the 
values range betwee n 0.8 8 and 0.99, depending on the details of the model, but 
remains below unity |l25| . We shall see (see Tables of numerical results) that axial 
coupling constant gA strongly depends (always grows with increasing dimension of 
representation) on the representation employed. In self-consistent quantum formal- 
ism in higher representations (j > 1/2) gA is always above unity even for small 
(^0.5 MeV) parameter f v values. 



3.4. Radii. In addition to electric mean square radius ( III. 3.12 ) the following 
three mean square radii are used: electric isovector, isoscalar magnetic, isovector 
magnetic. Below we give definitions for all of them following Ref. S|. 
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-0.2 



Figure HI. 3. Neutron charge density distribution. Short-dashed 
line denotes classical chiral angle distribution, with parameters 
e = 5.45; = 64.5 MeV from Ref. Quantum chiral an- 
gle distributions for various j values are plotted with solid lines. 
Long -dash ed line denotes semi-empirical charge density distribu- 
tion ||123||. 



For nucleons the isovector (electric) charge mean square radius becomes 

( r EJ=l) = -^p ( r E,I=l) 



1 J drr 4 sin 2 F^l + F 
e2 /' jdff 2 sin 2 F(l + F' 2 



12 i sin 2 _F 



i 2 F 



( r M,/=o) - -^2~f2 ( r MJ=o) 



The isoscalar magnetic mean square radius is expressed as 

1 

e 2 / 2 

12/ dff 4 F' sin 2 F 
~ ^Tf ^ fd?r 2 F< sin 2 F' 

and the isovector magnetic mean square radius takes the form 

( r M,7=l) = -^2~p ( r M,/=l) 

1 /d rT 4 S in 2 F(l + ^ + ^ + 4|^f(|-^ 
" e2 / 2 J dff 2 sin 2 F(l + F'2 + sis!£ + 4|^(| - 
All these radii can be measured experimentally. 
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(III.3.20) 



(III.3.21) 
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Table III . 1 . The predicted static baryon observables as obtained 
with the quantum Skyrme model for representations of different 
dimension. The first column (ANW) are the predictions for the 
classical Skyrme model given in Ref. |Q. Classical chiral angle 
function (taken from ANW column ) ha s been used for evaluation 
of integrals. The empirical results [127] are listed in the last col- 
umn. 





ANW 


j = l/2 


j = l 


j = 3/2 


j = 2 


j = 5/2 


Expt. 


m N 


Input 


Input 


Input 


Input 


Input 


Input 


939 MeV 




Input 


Input 


Input 


Input 


Input 


Input 


1232 MeV 


U 


64.5 


72.1 


76.4 


82.2 


89.4 


98.0 


93 MeV 


e 


5.45 


5.23 


5.15 


5.03 


4.89 


4.74 




ro 


0.59 


0.55 


0.53 


0.51 


0.48 


0.45 


0.72 fm 


M P 


1.87 


1.90 


1.84 


1.78 


1.71 


1.64 


2.79 




-1.31 


-1.42 


-1.40 


-1.37 


-1.35 


-1.33 


- 1.91 


9A 


0.61 


0.62 


0.65 


0.68 


0.73 


0.78 


1.26 


MA++ 




3.70 


3.58 


3.44 


3.29 


3.15 


3.7-7.50 


MA + 




1.71 


1.64 


1.55 


1.46 


1.37 


? 


Mao 




-0.28 


-0.31 


-0.34 


-0.38 


-0.42 


? 


Ma- 




-2.27 


-2.25 


-2.23 


-2.21 


-2.20 


? 



4. The static properties of the nucleon and A33-resonance in classical 

chiral angle approach 

The normalized state vector with equal fixed spin and isospin i are to be iden- 
tified with nucleons [l = ^) and A33 resonances {I — |). Before we proceed with 
the self-consistent quantum case it is wor th to examine quantum formulae with 
classical chiral angle solution (see Fig. [II. 2| ). We restrict this numerical analysis to 
irreducible representations only. Clearly, numerical values for reducible represen- 
tations lie in the range between minimal and maximal j values of which reducible 
representation is constructed. 

Classical chiral angle is obtained by solving the classical equation of motion 
( [I.1.30D (with appropriate boundary conditions) that is given by the requirement 
that the classical mass ( III.2.4l| ) be stationary. Asymptotic behaviour when f — > 00 
of classical solution can be easily found from asymptotic equation 



F"(r) + ?F'(f) - ^F(f) = 0. 



(III.4.1) 



Physical solution (satisfying F(f — > 00) = 0) of ( III.4.1 ) is 



k 

F(r-+ 00) = — 



k = const, 



(III.4.2) 



where k is determined by derivative continuity requirement and equation ( II.1.30 ) 
solution value at f = 0; ^(0) = it. Classical equations of motion can be solved, for 
example, in the following way. 



Recent measurements [128 



obtain value A 1 / \++ = 4.52. 
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Table III. 2. The predicted static baryon observables for rep- 
resentations of different dimension with fixed empirical values for 
the isoscalar radius and the axial coupling constant. 





j = 1/2 


j = 1 


j = 2 


j = 3 


j = 4 


j = 9/2 


3 = 7 


m N 


1434 


1402 


1300 


1147 


942 


821 


23 




1552 


1520 


1418 


1265 


1060 


939 


141 


u 


76.2 


76.2 


76.2 


76.2 


76.2 


76.2 


76.2 


e 


3.81 


3.81 


3.81 


3.81 


3.81 


3.81 


3.81 


ro 


Input 


Input 


Input 


Input 


Input 


Input 


Input 




4.18 


4.15 


4.07 


3.96 


3.80 


3.70 


3.08 


[In 


-3.85 


-3.83 


-3.75 


-3.63 


-3.47 


-3.38 


-2.76 


9A 


Input 


Input 


Input 


Input 


Input 


Input 


Input 




7.72 


7.67 


7.53 


7.32 


7.03 


6.86 


5.75 




2.90 


2.89 


2.84 


2.77 


2.67 


2.61 


2.24 


Ma° 


-1.92 


-1.90 


-1.85 


-1.78 


-1.69 


-1.63 


-1.26 


Ma- 


-6.73 


-6.69 


-6.55 


-6.33 


-6.05 


-5.88 


-4.77 



1. Fix a merge point r (it is assumed that asymptotic ( III.4.2|) giv es a good 
approximation to the solution at that point) and from ( [II. 4.2 ) calculate 
function F(fn) and its derivative F'(fo) values at that point. 

2. Choose arbitrary k and start standard differential equation iteration proce- 
dure until pointful r = is reached. 

3. If F(0) ^ 7r adjust k value, recalculate F(fo) and F'(fo) values and start 
differential equation iteration procedure again until F(0) = n to the required 
precision is achieved^. 

When classical solution is obtained, the c orrespon din g values for the Lagran gian 
parameters can be extracted from equations ( [II.4.3|) or ( III.3.17| ) and ( [II. 3.12 ), or 
from any of two combination s of the m. As in Ref. |8|, we determine the two 
parameters in the Lagrangian ( EII.2.3S| ) so that nucleon and A33-resonance masses 
t ake the ir empirical values. The expressions for these masses are then given by 
( HT.2.5SD 



m N = —M(F) 

e 

m A = —M(F) 
e 



J /,AM 3 (F) 



e 3 UAMj(F) 



2 • 45(F) ' 
15e 3 JV 



(III.4.3a) 



(III.4.3b) 



2 • 45(F) 

In the evaluation of these two masses numeric ally we employ the classical chiral 
angle F(r), t he sh ape of which is shown in Fig. [II. % 

In Table III. 1 we include the predicted values for static nucleon properties Q, 
as well as the original predictions obtained in Ref. [@] for the classical Skyrme model. 



28 Actually at f = equation ( 1I.1.3C ) becomes indefinite due to spherical coordinate system 
singularity at the origin. But we can solve the equation until some small value e. The function 
value at f = then can be calculated using Taylor-series expansion near the f = e point. From 
practical point of view it is convenient to solve a system of two first-order differential equations 
instead of one of the second order. 

29 There are, of course, much more efficient algorith ms fo r boundary value problems, for 
example, free and well known collocation software package p2q] (FORTRAN). 
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For larger values of j the quantum corrections become increasingly important. The 
key qualitative feature is that the quantum mass correction AMj(F) is negative, 
with a magnitude that grows with the dimension of the representation. Therefore it 
becomes possible to reproduce the empirical nucleon and A 33 -resonance mass values 
with increasingly realistic values of the pion decay constant /„■ as the dimension 
increases. This reaches its empirical value for a representation j — 5/2 of dimension 
6. There is an accompanying — if less significant — improvement of the numerical 
value for the axial coupling constant qa- 

In the case of the isoscalar radius r (Tables [[II. 1 and [II. 2) of the baryon, 
there is, however, no reduction of the difference between the predicted and the 
empirical value with increasing dimension of the representation. The same is true 
for the magnetic moments. The predicted value for the ratio of the proton and 
neutron magnetic moments actu ally d eteriorates slowly with increasing dimension 
of the representation. In Table [II.2 we show the representation dependence of 
the observables for the case where the parameters are determined by matching 
the empirical values of is oscalar r adius i~ej=o a nd the axial coupling constant gA- 
Because the expressions ( HI.3.17 ) and ( HI, 3.12] ) are independent of the dimension 
of the representation, parameters f n and e remain constant. The best agreement 
with the empirical values of the static properties of baryons in this case are obtained 
with the representation j = 4. 

We see that when the Skyrme model is treated consistently quantum mechan- 
ically ab initio the dimension of the representation becomes a significant addi- 
tional model parameter H|. By choosing two parameters of the model so as 
to match the empirical nucleon and A33-resonance masses it becomes possible to 
obtain a value for the pion decay constant, which is very close to the empirical 
value (89.4 MeV vs. 93 MeV). There was, however, no comparable gain in quality 
of the predictions for the baryon magnetic moments, which deteriorated slowly with 
increasing dimension of the representation. The value of axial coupling constant 
does on the other hand improve with increasing dimension, but stay below unity 
for representations of reasonably low dimension. 

When parameters of the model are chosen to match isoscalar radius and axial 
coupling constant, the parameters /„• and e are constants and nucleon and Ass- 
resonance observables depend only on the dimension of the representation. In con- 
trast to the first match, the magnetic moments improved with increasing dimension 
of the representation. Note that the treatment used here for the quantum skyrmion 
breaks down when the dimension of the representation grows so large that the neg- 
ative quantum mass correction AM,- becomes of the same ord er of magnitude as 
or larger than the classical skyrmion mass. As s hown in Table III.2 the numerical 
value of the quantum mass correction AMj(F) ( HI.2.42 ) is of the order of 55 MeV 
in the fundamental representation, but it rapidly increases in magnitude as the 
dimension of the representation grows. For a representation of dimension 9 (j = 4) 
it is large enough to cancel the ~ 500 MeV overprediction of the nucleon mass that 
is obtained when the empirical value for the pion decay constant is employed in the 
classical Skyrme model. For j — 15/2 the quantum mass correction exceeds the 
skyrmion mass and baryon masses become negative. After these comments we skip 
directly to quantum self-consistent treatment. 
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5. Self-consistent quantum formalism 



Minimization of the expression for the classical mass M(F) ( QI.2,41 ) leads 



to the conventional differential equation for the chiral angle ( II. 1.30 ) according to 
which F(f) falls as 1/f 2 at large distances. The behaviour is typical for long- 
range interaction and, therefore, implies zero pion mass. This is inconsistent with 
strong interaction properties, which is known to be short-range and according to 
Yukawa [ 128| | imply 

, . int. const ~ _ TTLC , . 

7r(r->oo = e~ mr ; rh = — , (III.5.1) 

r n 

pion field 7r(r) fall and, therefore, finite pion mass. T.H.R. Skyrme in his 1962 
work m wrote "These mesons have z ero ma ss, ultimately on account of the full ro- 
tational symmetry of the Lagrangian ( Il.l.g ). This symmetry is, however, destroyed 
by the boundary condition U(oo) = 1, and we believe that the mass may arise as 
a self consistent quantal effect." Thus, T.H.R. Skyrme noted that chiral group 
SU(2)l ® SU(2)r is spontaneously broken by choosing vacuum state U(oo) = 1 
to the subgroup diag(SU(2)L ® SU(2)r) ~ SU(2)j. According to the recent point 
of view just the spontaneous breakdown of chiral symmetry implies massless pions. 
In other words, pion is massless provided that the axial current is perfectly con- 
served. Because a pion is massive, we expect that axial symmetry should be an 
approximate symmetry and that the axial current should be only approximately 
(partially) conserved. 

In the semiclassi cal appr oach the quantum mass term AMe^ is absent from 
the mass expression ( III.2.52 ). The absence of negative AMs j correction has the 
consequence that va riati on of the truncated quantum mass expression yields no 
stable solution [34, 130 1 . The semiclassical approach describes the skyrmion as 
a "rotating" rigid body wi th fixed F(r) [^). In contrast the variation of the full 
energy expression ( III.2.52 ) that is obtained in the consistent canonical quantization 
procedure [48] in collective coordinates approach gi ves stab le solutions. 

Minimization of the quantum mass expression ( III. 2.52 ) leads to the following 
integro-differential equation for the chiral angle F(f) in the dimensionless variable 
f: 



F" 



+F'- 



-2r~ 2 



4sin 2 F 



-2sin2F 



15a 2 

g 4 ~2 gm 2JT 



■ (80aAM Ej + 20e(£ + 1) + 4di - 8d 2 sin 2 F 



15a 



7,1 



40dAAf Sj + 10£(£ + 1) + 2d x - 8d 2 sin 2 F 



F' -4f 



e rsm b 



+ sin2F 2 + 2 



15a 2 
sin 2 F 



(l60aAAf Sj + m{l + 1) + 8di - 16d 2 sin^F 

e 4 / 

^ ((40dAM s 



15a 2 



4d 3 sin 2 F 



1))(? 



lir 2 sin 2 ^ 



2 sin 2 F) 



15f 2 



0. 



(III.5.2) 

with the usual boundary conditions F(0) = ir and F(oo) — 0. The state dependence 
of this equation is a dir ect consequence of the fact that quan tization preceded 

variation (cf. Ref. | 131 |). Contrary to the classical equation ( II. 1 .3C| ) quantum 

chiral angle equation (111.5.2) depends on the parameter e. This, however, is not 
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an unusual result: quantization of breathing modes [3^, 37 1 also leads to parameter 
e dependent equations. 

At large distances this equation reduces to the asymptotic form 

r 2 F" + 2rF' - (2 + mlf 2 )F = 0, (III.5.3) 
where the quantity m\ is defined as 



ml 



3a(F) 



8AM S .(F) 



21(1 + 1) 



and the corresponding asymptotic solution takes the form 



F(f) 



l 



const. 



(III.5.4) 



(III.5.5) 



The requ ireme nt of sta bility of the quantum skyrmion is that the integrals ( |III.2.8| ) , 
( |TII2.41| ) and (|III.2.42|) converge. This requirement is satisfied only if m\ > 0. For 
that the presence of the negative quantum correction AM^ j (F) is necessary. The 
absence of this term leads to the instability of the skyrmion in the semiclassical 
approach [3_4|. Note that in the qua ntum t reatment the chiral angle possesses fi8| 
the asymptotic Yukawa behaviour ( |[II.5.5 ). The positive quantity m„ = ef^m^, 
therefore, can be interpreted as an effective mass for the pion field. 

It is known that in classical Skyrme model 3-div ergence of axial current density 
component £/ a (x) gives the differential equation ( II.1.3C ) for chiral angle |35|. 
"Therefore, in this chiral-symmetric theory the conservation of axial- vector current 
is, in fact, the equation of motion, i.e. the Euler-Lagrange equation for the rotated 
pion field tt," | pl| . Let's assume that in quantum theory axial symmetry is broken 
and the matrix element of the divergence of axial current is proportional to pion 
mass (PCAC relation^) 



£ 



£ 



f>lF(r)- 



(III.5.6) 



As a consequence we get the same asymptotic equation ( [II. 5.3] ), when r — > oo. 
This result supports the interpretation of m n as the effective pion mass. 

Note also that the finite pion m ass in the Skyrme model is usually introduced 
by adding an explicit term 48, fl32[ 

■S^ = Tr i U + Ut - 2 }' ( m - 5 - 7 ) 



in the Lagrangian density (11.1.6), and, therefore, the Lagrangian ( H.1.6| ) chiral 
symmetry becomes explicitly broken even in classical limit. 

Positivity of the pion mass ( III.5.4| ) can obviously only be achieved for states 
with sufficiently small values of spin t. This implies that the spectrum of states with 
equal spin and isospin will necessarily terminate at some finite value of the spin 
quantum number £. The termination point depends on parameters e and j values. 
When spin I value increases, the upper e value for which stable soliton solution 
exist always decreases (for the same repre sentation j). As the negative quantum 
mass correction AMj, in the expression ( QI.2,42 ) grows in magnitude with the 



To make explicit com parison with ( III. 3..^ ) just account that F(r)x a = j^-K a (see £[LLE )). 
Generally the comparison of (HI. 3. 3) and ( HI.5.6[ ) is complicated by the fact that in ( III.5 6| ) axial 
vector current i s sa ndwiched between two hadronic states. We will not consider this question here, 
but refer to j^, |3^ | and references therein. 
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dimension of the representation, it is always possible to find a representation in 
which the nucleon and the A33-resonance are the only stable particles, as required 
by experiment. 

Integro-differential equation ( III.5.2 ) can be attacked in the following way. 

1. Using chiral classical angle Fig. II. 2| and an y pair of empirical ba ryon obs erv- 
ables, for example, nucleon mass ( [II. 4.3a ) and isoscalar radius ( III.3.12| ) we 
fit two model para meters f n and e an d calcul ate all req ui red integ ra ls in th e 
quantum equation ( |III.5.2| ), namel y, fllll.2.41] ), ( [III.2.8D , ( |III.2.42| ), ( [111.5.41 ). 

2. Using known asymptotic solution ( HI.5.5 ) (and its derivative) one can adopt 
simple differential equation solution procedure described in Sec. |] and find 
the first quantum solution F\(f) and the constant k\. 

3. This quantum solution F\{ f) now can be subsequently used to recalculate 
U,e and integrals ( [HI.2 4l| ), ( [HI.2.8D , ( |III.2.42| ), ( |III.5.4| ). Then again pro- 
cedure described in item |2j can be used to get the solution ^(f) and constant 
k 2 . 

4. This procedure can be iterated until convergent solution and parameters 
/tt, e as well as stable values of M(F), AMj^ (F), a(F), m,. K (F) are obtained. 
The self-consistent set then can be used to calculate numerous phenomeno- 
logically interesting quantities. 

Quant um chir al angle solutions wit h mode l parameters dete rmine d from nucleon 
mass ( III.4.3a|) and isoscalar radius ( HI.3.12 ) are shown in Fig. III. 4 . If one succeeds 
in initial guessQ, then it usually takes only 10-15 iterations to get 5-6 fixed digits 
in all integrals and f„,e values. For nucleon solutions with spin £ = 1/2 exist (in 
the fundamental representation j — 1/2) when < e < 7.5. The largest value of 
e for which stable solutions are obtained decreases with increasing dimensionality 
of the representation, and there are no restrictions to the existen ce of s olution 
from representation j employed. This can b e seen from equation ( [II. 5.2 ) itself. 
I ndeed, the formal substitution e — > into ( [II. 5.2 ) yields the clas sical equation 
( [1.1.30 ), which does have a solution^. For illustra tion (see Fig. III.5 ) we also 
provide nucleo n ma ss densities (Hamiltonian densities ( [II. 2.52 )) for quantum chiral 
angles in Fig. [11.4 and detailed contributions of classical soliton masses, rotation 
ener gies a nd quantum mass corrections for two of them (j = \ and j 
Fig. " 



[II.6 



A remark is given t o the a model without the Skyrme term (this corresponds 
to the limit e — ► oo in ( H.1,6] )). Although the quantum correction AM^ plays a 
role in stabilization of soliton, one cannot obtain "quantum" rotating chiral soli- 
ton without the Skyrme term as stabilizer | |133|| . Applying higher representations 
also cannot help in obtaining this stability |134| | . This is consistent with Derrick 
theorem [57]. For A33-resonances {£ — 3/2) there are no stable solutions in the 
fundamental (j = 1/2) representation, nor in the representation j = 1. In the rep- 
resentations j = 3/2 and j = 2 there are only stable soliton solutions for baryons 
with spin I = 1/2 and I = 3/2. A dimension with j = 5/2 allows solitons with 
£ = 1/2, 3/2 and 5/2, and, therefore, appears to be empirically contraindicated. 
The numerical results for nuc leon s with fixed empirical values for isoscalar radius 
and mass are shown in Table A^2 in Appendix [a|. For the irreducible representa- 
tion j = 1 the proton magnetic moment calculated in this way is within 10% of 

31 Solution of classical equa tion I H.1.3o| ) exist or not for all e, f n values. This is not the case 
for integro-differential equation ( HI.5.2| f! 

32 Note that e — > does not lead to any meaningful limit due to dimensionless variable 
f = ef^r being used. This, however, does not affect the solution existence argument for arbitrary 
small e, because only formal similarity of equations is important. 
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Figure III. 5. Family of nucleon mass densities in representations 
j = 2 j 1 ffi 2 ® 2' ^' §» § ( m °del parameters being calc ulated fr om 
empirical nucleon mass ( III, 4.3a ) and isoscalar radius ( III.3. 12| )) . 
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Figure III. 6. Detailed contributions of classical soliton masses, 
rotation energies and quantum mass corrections for nucleon in rep- 
resentations j = \ and j=l (mo del para meters being calculated 
from em pirical nucleon mass ( III. 4.3a ) and isoscalar radius 

mi)). 



the empirical value. The calculated values of both the neutron magnetic moment 
and the axial coupling constant agree with the corresponding empirical values to 
within 1%. Numerical r esult s for fixed isoscalar radius and axial coupling constant 
are represented in Table |A.3| (Appendix [a]) . For representation j = 1 nucleon mass 
and proton magnetic moment agree wit h ex perimental valu es within 10%, neutron 
magnetic moment within 2%. In Table |A.l| and Table |A.^ we rep rese nt numerical 
results obtain ed b y fixing — 138 MeV and nucleon mass (Table AA) or isoscalar 
radius (Table AA). In these approaches agreement with experimental data is worse. 
We note, however, that pion field mass used here (138 MeV) approximately equals 
to particles 7ro,7r_,7r + mass. The effective strong interaction field fall can be in- 
directly estimatecp^| from experiments. These measurements corresp ond to about 
150 MeV effective boson mass, when extracted from Yukawa formula ( HI. 5.1 ). 

For parameters e, f n fixed from nucleon data, empirical values for A33-resonance 
observables for different representations can be obtained. Results for e = 4.15, = 
58.5 MeV (determined by fit to nucleon m N = 939 MeV and (r| /=0 ) 1/2 = 0.72 fm 



and representation j = 1) are demonstrated in Table [B.l| in Appendix |Bj. Re- 
sults for the same fit o f experimental qua ntitie s, but in different repres entations 
are shown in Table [|]| (j = \) and Table IjO (j = 1 | © h)- Tables |B^, [R7, 



B.8 represent numerical results with fixed values for e and f^, which are calculated 



from nucleon observables gA and 



J = l 



EJ=0/ 



a/2 



in representations j 



j = 1 and 



s, respectively. Tables B.4 and B.5 show similar results for fixed e, 



Private communications of E. Norvaisas and D.O. Riska. 



62 



III. QUANTUM SKYRME MODEL 



values, extracted from nucleon observables and {r 2 E J=0 ) 1 ^ 2 in representations 
j = 1 (Table |D|) and j = 1 ± § (Table pPf) . 



It is also worth to mention the numerical results for reducible representations 
1 <8> \ <S> ^ for nucleon and | <£> 1 <8> 5 for A33-resonance. The representation § <£> 1 <8> \ 
is the first one for which the stable soliton solution is obtained for I = 3/2 state. 
The numerical results in these reducible representations are comparable — if not 
better — to those of representations j = 1 (for £ = 1/2) and j = 3/2 (for £ = 3/2), 
respectively. Noting that the SU(3) isomultiplets — octet (1, 1) and decuplet (2, 0) 
— (to which nucleon and A33-resonance belong, respectively) split, when restricted 
to SU(2) subgroup, into 1 <g> | <g> \ <g> and | ® 1 ® | ® 0, respectively^, the following 
explanation is not improbable: it just can indicate that SU(2) Skyrme model in 
some way (via representations employed) reveals relevant SU(3) symmetry of strong 
interactions. If this is indeed the case, then in SU(3) Skyrme model using higher 
representations ((1, 1) and (2,0)) one can hope to describe the entire isomultiplets 
of baryons more successfully. Investigation of SU(3) Skyrme model is required to 
confirm or reject these considerations. 

6. Remarks on persisting problems 

Despite quite successful generalization of quantum Skyrme model to general 
representation of SU(2) group a few problems still persist. First of all we see rather 
big gap between calculated (/ T ~ 60 MeV, e ~ 4) and extracted from experiment 
(./tt = 93 MeV, e ~ 7.4) numerical values of model parameters. The problem is 
inherited from semiclassical calculations |§. The difference, however, becomes in- 
creasingly important, because all other measurable quantities fit to their experimen- 
tal values much better (~ 10%, see Appendix ^). Thus the question again arises to 
what extent do these errors of the Skyrme model reflect use of the approximation, 
and to what extent do they reflect the fact that the Skyrme model is crude approx- 
imation to meson physics. The question has been considered by E. Witten who 
suspected that in the case of semiclassical approximation "the error mainly has the 
latter origin. . . If the proper meson theory were known, the semiclassical approx- 
imation would be equivalent to the 1/N C expansion, and I personally suspect the 
error would be much less than 30%". Our calculations, if we exclude f v ,e values, 
which can be considered just as model parameters without referring to its physical 
content, are in much better agreement with experimental data. However, it should 
be noted that we obtain these values due to an additional implicit discrete param- 
eter — the group representation. The quantization procedure also differs from the 
semiclassical approach |^| considered by E. Witten J^|. A problem is also present 
with the interpretation of representation dependence itself. Despite the quantiza- 
tion method used here is free of ambiguity and is in agreement with principles of 
modern quantization of constrained systems, it undoubtedly shows representation 
dependence of measurable quantities^ Thus physical interpretation is welcome 
and we hope that quantum investigation of SU(3) Skyrme model can enlighten this 
problempt 

34 Representation j = does n ot influ ence the results. Indeed, j = representation of SU(2) 



is trivial and, thus, the Lagrangian (II.1.6|) is identical zero: S£ = 



35 The problem we suspect is not the problem of the Skyrme model only but rather of the 
quantization procedure itself. As a consequence, we expect representation dependence in other 
models defined on group manifold. 

36 We hope that the situation is somehow similar to that of W. Heisenberg and M. Gell-Mann 
about assumptions of symmetry of strong interaction, briefly mentioned in the Preface. 



Concluding statements 



The following statements represent the main results of the work in the order of 
decreasing importance. 

• Each of SU(2) representation j yields the different quantum Lagrangian 
density. As a consequence, theoretical observables depend on representation 
j which can be treated as a new phenomenological parameter. 

• Quantum chiral solitons exist and po ssess a symptotic behaviour consistent 
with the massive Yukawa field fall ( III. 5.1 ). The asy mptotic shape and 
PC AC relation leads to the correct asympto tic equa tion ( QI.5.3 ) coinciding 
with contribution of explicitly broken term ( [II. 5.7 ). This encourages us to 
suggest that the integral m w should be interpreted as an effective pion mass. 

• A nucleon and A33-resonance are the only stable states for irreducible rep- 
resentations j = | and j = 2. Unphysical tower of states £ sp ; n = Asospim 
which is artifact of the classical and semiclassical Skyrme model is, therefore, 
terminated by choosing the appropriate SU(2) representations. 

• Higher spin (£ > 1/2) quantum state s are not "spherically symmetric". The 
Hamiltonian ( III.2.53 ) (Lagrangian ( III.2.4C| )) density function depends on 
the polar angle Nucleon states are "spherically symmetric" in various 
representations of j, as required. 

• Each of the spin-isospin state yields the different range of realizable values 
of the parameter e. A stable quantum self-consistent solution exist: 

— for spin t—\ states in all SU(2) representations, 

— for I = § spin states starting at least from reducible representation 

3 m 1 rr, 1 

2 w w 2 " 

• A very good agreement with experimental data is obtained for axial coupling 
constant gA in higher representations of SU(2), the problem being previously 
unsolved by using various extensions of the model in the fundamental rep- 
resentation of SU(2). 

• The basic quantum Skyrme model provides considerable improvements in 
nucleon magnetic momenta and, especially, in neutron charge density distri- 
bution. 

The text of PhD thesis can be found at the Website of Institute of Theoretical 
Physics and Astronomy 
http:/ /www. itpa.lt/baryorj 



CONCLUDING STATEMENTS 



APPENDIX A 



Nucleon observables in different representations 



Table A.l. The predicted static baryon observables for different 
representations with fixed empirical values for the effective pion 
mass m w = 138 MeV and nucleon mass 939 MeV. 



j 


1/2 


1 


3/2 




Expt. 


m N 


Input 


Input 


Input 


Input 


939 MeV 


u 


68.4 


54.9 


49.6 


57.2 


93 MeV 


e 


4.97 


3.96 


3.52 


4.15 






1.63 


3.20 


4.44 


2.80 


2.79 




-1.06 


-2.61 


-3.85 


-2.21 


-1.91 


9A 


0.89 


1.43 


1.77 


1.30 


1.26 


m w 


Input 


Input 


Input 


Input 


138 MeV 


\j ( r E,I=o) 


0.54 


0.83 


1.00 


0.77 


0.72 fm 



Table A. 2. The predicted static nucleon observables in different 
representations with fixed empirical values for the isoscalar radius 
( r l./=o) 1/2 = °- 72 fm and nucleon mass 939 MeV. 



j 


1/2 


1 


3/2 


5/2 




Expt. 


m N 


Input 


Input 


Input 


Input 


Input 


939 MeV 


u 


59.8 


58.5 


57.7 


56.6 


58.8 


93 MeV 


e 


4.46 


4.15 


3.86 


3.41 


4.24 






2.60 


2.52 


2.51 


2.52 


2.53 


2.79 


Mn 


-2.01 


-1.93 


-1.97 


-2.05 


-1.93 


-1.91 


9A 


1.20 


1.25 


1.33 


1.52 


1.23 


1.26 




79.5 


180. 


248. 


336. 


155. 


138 MeV 


\J ( r E,I=o) 


Input 


Input 


Input 


Input 


Input 


0.72 fm 


y/(r%,i=i) 


1.33 


1.03 


0.97 


0.93 


1.07 


0.88 fm 


\J ( r M,7=o) 


1.05 


1.01 


1.00 


1.00 


1.01 


0.81 fm 


\J ( r M,7=l) 


1.32 


1.03 


0.97 


0.93 


1.07 


0.80 fm 
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A. NUCLEON OBSERVABLES IN DIFFERENT REPRESENTATIONS 



Table A. 3. The predicted static baryon observables for different 
representations with fixed empirical values for the isoscalar radius 
(r% /=0 ) 1 ^ 2 = 0.72 fm and axial coupling constant gA — 1-26 . 



J 


1/2 


1 


3/2 


5/2 




Expt. 


m N 


986. 


948. 


882. 


694. 


963. 


939 MeV 


U 


61.4 


58.9 


55.7 


48.1 


59.7 


93 MeV 


e 


4.37 


4.13 


3.92 


3.56 


4.20 




Mp 


2.70 


2.53 


2.42 


2.05 


2.57 


2.79 




-2.14 


-1.95 


-1.84 


-1.51 


-1.99 


-1.91 


9A 


Input 


Input 


Input 


Input 


Input 


1.26 




75.8 


179. 


259. 


386. 


152. 


138 MeV 


\J ( r E,I=o) 


Input 


Input 


Input 


Input 


Input 


0.72 fm 


l/(r 2 EJ=1 ) 


1.36 


1.03 


0.96 


0.92 


1.08 


0.88 fm 


\] ( r MJ=o) 


1.05 


1.01 


1.00 


1.01 


1.01 


0.81 fm 


\J ( r M,7=l) 


1.35 


1.03 


0.96 


0.92 


1.08 


0.80 fm 



Table A. 4. The predicted static baryon observables for different 
representations with fixed empirical values for the isoscalar radius 
i r E /=o) 1 ^ 2 = 0-72 f m an d effective pion mass — 138 MeV. 



J 


1/2 


1 


3/2 


5/2 




Expt. 


m N 


486. 


1264. 


1872. 


2992. 


1069. 


939 MeV 


U 


41.7 


68.8 


84.0 


106.3 


63.1 


93 MeV 


e 


5.56 


3.76 


3.12 


2.48 


4.05 




Mp 


1.93 


3.13 


4.35 


6.72 


2.77 


2.79 




-0.86 


-2.68 


-4.04 


-6.52 


-2.24 


-1.91 


9A 


0.60 


1.65 


2.46 


3.94 


1.39 


1.26 


m T 


Input 


Input 


Input 


Input 


Input 


138 MeV 


\J ( r E,I=o) 


Input 


Input 


Input 


Input 


Input 


0.72 fm 


l/(r 2 EJ=1 ) 


1.07 


1.12 


1.12 


1.12 


1.11 


0.88 fm 


\J ( r MJ=o) 


1.01 


1.02 


1.02 


1.02 


1.02 


0.81 fm 


\J ( r M,7=l) 


1.05 


1.12 


1.12 


1.13 


1.11 


0.80 fm 



APPENDIX B 



A33-resonance observables in different 
representations 



Table B.l. The predicted static A33-resonance observables in dif- 
ferent representations with fixed values for the parameters e = 4.15 
and fn — 58.5 MeV (from nucleon observables tojv = 939 MeV, 



E,I=0 



< x < 2 = 0.72 fm, representation j = 1 in Table A.2). 



j 


§©i©^ 


a 

2 


2 


Expt. 




1055. 


1029. 


910. 


1232 MeV 


MA++ 


7.38 


6.40 


4.20 


3.7-7.5 




3.02 


2.73 


2.01 


? 




-1.33 


-0.94 


-0.19 


? 


Ma- 


-5.69 


-4.61 


-2.38 


? 


\J ( r E,I=o) 


0.91 


0.87 


0.72 


? 



Table B.2. The predicted static A33-resonance observables for 
different representations with fixed empirical values for the e = 
4.46 and — 59.8 MeV (from nucleon observables tojv = 939 MeV, 
( r i,/=o) 1/2 = °- 72 fm > representation j = \ in Table [IJ). 



j 


3 m -i m 1 


2 


2 


Expt. 




1008. 


974. 


809. 


1232 MeV 


MA++ 


6.05 


5.15 


3.00 


3.7-7.5 


MA+ 


2.63 


2.36 


1.63 


? 




-0.80 


-0.43 


0.25 


? 


Ma- 


-4.23 


-3.22 


-1.12 


? 




104. 


172. 


438. 


138 MeV 


\/ ( r E,i=o) 


0.84 


0.79 


0.62 


? 
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Table B.3. The predicted static A33-resonance observables for 
different representations with fixed empirical values for the e = 
4.24 and f n = 58.8 MeV (from nucleon observables m N = 939 MeV, 
( r l./=o) 1/2 = °- 72 fm > representation j = 1 © \ © \ in Table [Ofr . 



j 




a 
2 


2 


Expt. 




1040. 


1012. 


881. 


1232 MeV 


MA++ 


6.96 


6.00 


3.82 


3.7-7.5 


MA+ 


2.90 


2.61 


1.89 


? 


^ A 


-1.16 


-0.78 


-0.05 


? 


Ma- 


-5.23 


-4.17 


-1.98 


? 




85.5 


141. 


338. 


138 MeV 


\J ( r E,I=o) 


0.89 


0.85 


0.69 


? 



Table B.4. The predicted static A 33 -resonance observables for 
different representations with fixed empirical values for the e = 
3.76 and /„■ = 68.8 MeV (from nucleon observables = 138 MeV, 
( r Ej=o) 1/2 = °- 72 f m > representation j = 1 in Table |A~J). 



j 




a 
2 


2 


Expt. 




1360. 


1338. 


1245. 


1232 MeV 


MA + + 


8.14 


7.24 


5.28 


3.7-7.5 


MA+ 


3.14 


2.86 


2.23 


? 


Ma° 


-1.87 


-1.52 


-0.83 


? 


Ma- 


-6.87 


-5.89 


-3.88 


? 




65.3 


106.5 


234. 


138 MeV 


^/ ( r E,I=o) 


0.84 


0.81 


0.72 


? 



Table B.5. The predicted static A 33 -resonance observables for 
different representations with fixed empirical values for the e = 
4.05 and f v = 63.1 MeV (from nucleon observables = 138 MeV, 
( r l./=o) 1/2 = °- 72 frn i representation j = 1 © \ © \ in Table |A~J). 



j 


§©iffl§ 


2 


2 


Expt. 




1162. 


1137. 


1020. 


1232 MeV 


MA++ 


7.25 


6.32 


4.27 


3.7-7.5 


MA+ 


2.92 


2.65 


1.97 


? 


Ma 


-1.40 


-1.03 


-0.32 


? 


Ma- 


-5.72 


-4.71 


-2.62 


? 




78.5 


129. 


297. 


138 MeV 


\/ ( r l,/=o) 


0.86 


0.82 


0.70 


? 
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Table B.6. The predicted static A33-resonance observables for 
different representations with fixed empirical values for the e = 
4.37 and /„• = 61.4 MeV (from nucleon observable s qa = 1-26, 
( r l./=o) 1/2 = °- 72 f m > representation j = \ in Table |a]|). 



j 




a 

2 


2 


Expt. 




1053. 


1021. 


865. 


1232 MeV 


MA++ 


6.18 


5.28 


3.19 


3.7-7.5 




2.64 


2.37 


1.67 


? 




-0.90 


-0.54 


0.14 


? 


Ma- 


-4.44 


-3.45 


-1.39 


? 




99.7 


165. 


409. 


138 MeV 


\J ( r E,I=o) 


0.83 


0.79 


0.62 


? 



Table B.7. The predicted static A 33 -resonance observables for 
different representations with fixed empirical values for the e = 
4.13 and f„ = 58.9 MeV (from nucleon observables qa = 1.26, 
( r Ej=o) 1/2 = 0- 72 f m > representation j = 1 in Table [Of). 



j 


fffiiffi^ 


a 

2 


2 


Expt. 




1064. 


1038. 


920. 


1232 MeV 


MA + + 


7.40 


6.42 


4.23 


3.7-7.5 


MA+ 


3.03 


2.73 


2.01 


? 


Ma° 


-1.35 


-0.96 


-0.21 


? 


Ma- 


-5.72 


-4.65 


-2.43 


? 




78.4 


129. 


302. 


138 MeV 


\J ( r E,I=o) 


0.91 


0.87 


0.72 


? 



Table B.8. The predicted static A 33 -resonance observables for 
different representations with fixed empirical values for the e = 
4.20 and = 59.7 MeV (from nucleon observables qa — 1-26, 
( r l./=o) 1/2 = °- 72 fm ' representation j = 1 © ± © ± in Table [Aj|. 



j 


3 ffil <c 1 

2 w ■*- W ^ 


a 

2 


2 


Expt. 




1062. 


1035. 


907. 


1232 MeV 


MA++ 


7.01 


6.06 


3.90 


3.7-7.5 


MA+ 


2.90 


2.61 


1.90 


? 


Ma 


-1.21 


-0.83 


-0.10 


? 


Ma- 


-5.32 


-4.27 


-2.10 


? 




84.2 


139. 


330. 


138 MeV 


\/ ( r B,I=o) 


0.88 


0.84 


0.69 


? 
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APPENDIX C 

Baryon densities for B > 1 configurations 



c 







. <1 * . 



Figure C.l. Surfaces of constant baryon density for the following 
Skyrme fields |^|: a) B = 2 torus b) 5 = 3 tetrahedron c) 5 = 4 
cube d) £> = 5 with Did symmetry e) B — 6 with symmetry f ) 
B = 7 dodecahedron g) B — 8 with symmetry h) B — 9 with 
tetrahedral symmetry i) B — 17 buckyball j) £? = 5 octahedron k) 
B = 11 icosahedron. 
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